Temperature and bath size in exact diagonalization 
dynamical mean field theory 
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Dynamical mean field theory (DMFT) combined with finite-temperature exact diagonalization 
is one of the methods to describe electronic properties of strongly correlated materials. Because 
of the rapid growth of the Hilbert space, the size of the finite bath used to represent the infinite 
lattice is severely limited. In view of the increasing interest in the effect of multi-orbital and multi- 
site Coulomb correlations in transition metal oxides, high-T c cuprates, iron-based pnictides, organic 
crystals, etc., it is appropriate to explore the range of temperatures and bath sizes in which exact 
diagonalization provides accurate results for various system properties. On the one hand, the bath 
must be large enough to achieve a sufficiently dense level spacing, so that useful spectral information 
can be derived, especially close to the Fermi-level. On the other hand, for an adequate projection 
of the lattice Green's function onto a finite bath, the choice of the temperature is crucial. The 
role of these two key ingredients in exact diagonalization DMFT is discussed for a wide variety of 
systems in order to establish the domain of applicability of this approach. Three criteria are used to 
illustrate the accuracy of the results: (i) the convergence of the self-energy with bath size, (ii) quality 
of the discretization of the bath Green's function, and (iii) comparisons with complementary results 
obtained via continuous-time quantum Monte Carlo DMFT. The materials comprise a variety of 
three-orbital and five-orbital systems, as well as single-band Hubbard models for two-dimensional 
triangular, square and honeycomb lattices, where non-local Coulomb correlations are important. 
The main conclusion from these examples is that a larger number of correlated orbitals or sites 
requires a smaller number of bath levels. Down to temperatures of 5 to 10 meV (for typical band 
widths W ~ 2 eV) two bath levels per correlated impurity orbital or site are usually adequate. 

PACS. 71.27+a Strongly correlated electron systems 



Contents 

1. Introduction 

2. ED DMFT Formalism 

2.1. Multi-Band Correlations 

2.2. Multi-Site Correlations 

2.3. Size of Hilbert Space 

2.4. Bath Discretization 

3. Results and Discussion 

3.1. Degenerate Three-Band Model 

3.2. Three-Band Materials 

3.3. Degenerate Five-Band Model 

3.4. Triangular Lattice 

3.5. Square Lattice 

3.6. Honeycomb Lattice 

4. Summary and Outlook 
Acknowledgments 
Appendix: Bath Discretization 

A.l. Degenerate Three-Band Model 
A. 2. Square Lattice 
References 



1. INTRODUCTION 

The study of strong Coulomb correlations is one of the 
central topics in condensed matter physics. A wide vari- 



ety of phenomena, such as high-T c superconductivity in 
the cuprates, paramagnetic metal insulator transitions, 
strong effective mass enhancements, and bad-metallic be- 
havior are intimately related to electronic correlations 
not included in ordinary band theory. During the past 
two decades, dynamical mean field theory [3-Q (DMFT) 
has proved to be a highly useful scheme for the descrip- 
tion of these kinds of phenomena in many materials. The 
virtue of DMFT is that intricate single-electron proper- 
ties, associated with complex lattice geometries and unit 
cells, and many-electron interactions in partially filled, 
nearly localized atomic orbitals are treated on the same 
footing. Since the balance between kinetic energy and 
Coulomb repulsion depends in subtle ways on important 
parameters, such as temperature, pressure, and chemical 
doping, a sufficiently detailed treatment of all key ingre- 
dients is essential. The conceptual advance of DMFT 
is that the problem of the correlated infinite lattice is 
mapped onto the problem of a correlated impurity im- 
mersed in a bath whose self-energy is determined via a 
self-consistent iterative procedure. The effective impu- 
rity may consist of a single atom or a small cluster of 
atoms. Coulomb interactions within this impurity are 
treated numerically exactly. Thus, dynamical fluctua- 
tions, which give rise to spectral weight transfer between 
low and high energies, are fully taken into account. Spa- 
tial fluctuations, on the other hand, are either neglected, 
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as in local or single-site DMFT, or included only at short 
distances, as in cluster extensions of DMFT 

The main computational task in DMFT lies in the so- 
lution of the quantum impurity problem for which tech- 
niques such as numerical renormalization group IH 14 1 
(NRG), quantum Monte Carlo (QMC), exact di- 
agonalization [16] (ED), and other more approximate 
schemes can be used. For multi-band materials and clus- 
ters involving several atoms, the most accurate and ver- 
satile approaches are QMC and ED. The focus in the 
present work is on ED at finite temperatures 17H21 1 . 



The key idea in ED is to discretize the continuum of 
the infinite lattice via a finite set of bath levels and to 
solve the many-body problem involving correlated impu- 
rity and finite bath exactly. Originally applied to the 
one-band Hubbard model 16j, it soon was realized that 



this approach is limited to a small number of bath levels 
because of the exponentially growing Hilbert space. The 
usefulness of ED DMFT for strongly correlated materials 
involving several orbitals or sites was therefore far from 
evident. 

The purpose of this work is to demonstrate that finite- 
temperature ED DMFT is in fact a highly useful and ac- 
curate many-body scheme that is applicable to systems 
including (so far) up to five correlated orbitals and clus- 
ters up to six correlated sites. Typically, this implies two 
to three bath levels per orbital or site. In certain cases, 
only one bath level per orbital or site was found to give 
satisfactory results. This conclusion seems surprising in 
view of the fact that, to obtain well converged results for 
the single-band Hubbard model within local DMFT, at 
least three to five bath levels must be taken into account 
An analogous five-orbital system would imply a to- 
tal number of 20 to 30 levels, a virtually insurmountable 
and impractical task. The main point of our study is 
that an increasing number of correlated orbitals or sites 
requires a decreasing number of bath levels, so that, for 
the systems discussed in this work, overall sizes of 10 to 
15 levels are adequate. Since computational times per 
iteration are of the order of minutes for 10 levels and a 
few hours for 15 levels, this implies that a large variety 
of strongly correlated materials of interest can be studied 
within ED DMFT. Moreover, this approach is applicable 
at large Coulomb energies and for full Hund exchange. 
Also, it is free of sign problems and statistical errors. We 
therefore conclude that finite-temperature ED DMFT is 
complementary to continuous-time QMC DMFT (22H24 1 . 
While the latter becomes computationally more demand- 
ing at low temperature because of statistical noise, the 
former gets more involved at high temperature due to 
the increasing number of excited states. Moreover, the 
low-temperature, low-energy properties are affected by 
finite-size limitations. 

As discussed in Section 2, at each iteration step ED 
DMFT implies two projections: from the infinite lat- 
tice to the finite cluster, and back from the cluster to 



the lattice. Mathematically, these projections assume 
that continuous and discrete real-energy spectra of the 
self-energy and Green's function exist which have nearly 
identical representations along the Matsubara axis. To 
demonstrate the accuracy of this approach three criteria 
are available: (i) the convergence of key quantities, such 
as the self-energy, with increasing bath size, (ii) the ac- 
curacy of the discretization of the lattice bath Green's 
function in terms of a cluster consisting of impurity and 
finite bath (the impurity may be formed by several sites) , 
and (iii) the consistency with analogous results derived 
within QMC DMFT. In the case of three-band systems, 
calculations for two, three, and four bath levels per or- 
bital have been performed. The comparisons for different 
bath sizes, the quality of the Green's function discretiza- 
tion, and the comparison with QMC results illustrate 
that not only integrated quantities, such as orbital oc- 
cupancies, but also dynamical features are usually well 
described with only two bath levels per orbital or site. 

The successful application of ED DMFT involves two 
main ingredients: (i) the size of the Hilbert space and (ii) 
the discretization of the lattice via a finite cluster. The 
size of the Hilbert space determines the level spacing of 
excited states. This spacing must be sufficiently small, so 
that the low-energy behavior of the self-energy can be re- 
solved. Of particular interest are correlation-induced ef- 
fective mass enhancements, deviations from Fermi-liquid 
properties related to finite low-energy scattering rates, 
and the opening of Mott gaps. The accuracy of the 
discretization, on the other hand, depends in a crucial 
manner on the temperature used along the imaginary 
Matsubara axis. Decreasing this temperature implies 
increasing sensitivity to low-energy features of the lat- 
tice Green's function which require baths of accordingly 
larger size. As the cluster Green's function is always 
gapped, it cannot describe metallic behavior at very low 
energies. The examples discussed in the subsequent sec- 
tions indicate that, to achieve accurate results with two 
bath levels per orbital or site, the Matsubara tempera- 
ture should not be much lower than T w 5, . . . , 10 meV 
(assuming typical band widths W « 2 eV). 

In the past, several groups have carried out multi-band 
or multi-site ED DMFT calculations at T = [H^. In 
these studies, the impurity Green's function is evaluated 
using the Lanczos method for the ground state while the 
DMFT iteration procedure and the discretization of the 
bath Green's function are performed at a fictitious finite 
temperature to obtain sufficiently smooth distributions 
along the Matsubara axis. In the present work we focus 
mainly on ED at finite temperature in order to demon- 
strate the consistency of this approach with analogous 
results derived within continuous-time QMC DMFT. As 
in QMC, the temperature for the Matsubara grid is then 
the same as for the interacting impurity Green's func- 
tion. At low temperatures this function can in principle 
be obtained by generalizing the Lanczos scheme to finite 
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T. This approach, however, is rapidly plagued by the 
loss of orthogonality of excited states [21| . Here we make 
use of the Arnoldi algorithm [33| for large sparse ma- 
trices, which is ideally suited to evaluate those excited 
states that are physically important at low temperatures 
[20| . Whereas finite-T ED with full matrix diagonaliza- 
tion is feasible only up to about 8 levels [3] (maximum 
Hamiltonian dimension N — 4900), the Arnoldi scheme 
at present allows the consideration of 12 to 15 levels 
(N fx 10 6 , . . . , 40 x 10 ) without significantly increased 



storage requirements 34|. Thus, materials involving up 



to five correlated orbitals or up to six correlated sites can 
now be investigated. 

As the ED DMFT iterative procedure is carried out 
at imaginary frequencies, an additional step is required 
to extrapolate the converged Green's function and self- 
energy to real u). In principle, these quantities can, of 
course, be evaluated directly at 10 + £7, where 7 denotes 
an artificial broadening. This aspect represents an impor- 
tant advantage compared to QMC DMFT, where real- 
energy spectra are usually obtained via the Maximum 
Entropy method (35|. On the other hand, the ED cluster 
properties are by definition discrete, so that the imagi- 
nary parts of the Green's function and self-energy consist 
of (broadened) <5-functions. These cluster spectra are ad- 
equate, for instance, if a metallic phase with many closely 
spaced peaks is compared to an insulating phase, which 
exhibits a clear gap. Mott transitions can therefore read- 
ily be identified via cluster spectra. To investigate finer 
features, such as the distinction between Fermi-liquid and 
non-Fermi-liquid behavior, it is useful (as in QMC) to 
analytically continue the lattice quantities from the Mat- 
subara axis to real uj. Although this problem is known 
to be ill-posed, the important low-energy properties are 
usually rather insensitive to details of the extrapolation 
method, while large uncertainties exist in the incoherent 
region farther away from the Fermi energy. To illustrate 
some of these issues, cluster and lattice spectra will be 
compared in Section 3 for several systems. 

Multi-orbital and multi-site Coulomb correlations are 
treated in this work within a common perspective since 
they share many computational aspects. The similar- 
ity becomes particularly evident by changing from the 
site basis in the case of cluster DMFT to a "molecular- 
orbital" or "plaquette" basis in which the onsite repulsion 
is converted into intra-orbital and inter-orbital Coulomb 
and exchange interactions. An interesting outcome of 
this parallel treatment is that characteristic physical phe- 
nomena, such as the transition from a Mott insulating 
phase at half-filling to non-Fermi-liquid and Fermi-liquid 
metallic phases at finite electron or hole doping, are found 
in multi-orbital materials as well as single-band models 
of high-T c cuprates J^. Thus, the existence of inter- 
orbital or inter-site interaction channels generates fasci- 
nating new physics not present in a single-band, single- 
site picture. 



The outline of this paper is as follows. In Section 2 
we review the finite-temperature ED DMFT approach 
for the evaluation of multi-band and multi-site Coulomb 
correlations. We also discuss the role of the two main 
ingredients in this approach, namely, the overall size of 
the Hilbert space of the finite cluster used to represent 
the infinite lattice, and the discretization of the bath 
Green's function. Section 3 presents the results for a 
variety of applications: (1) the degenerate three-band 
model, (2) several realistic three-band materials, such as 
Ca 2 - x Sr :c Ru04 ) LaTi0 3 , V 2 3 , and Na x Co0 2 , (3) the 
degenerate five-band model and its connection to iron- 
based pnictides and chalcogenides, the single-band Hub- 
bard model for (4) the triangular lattice, (5) the square 
lattice, and (6) the honeycomb lattice. We have also ap- 
plied finite-T ED within inhomogeneous DMFT [HI to 
study both single-site and non-local correlations in vari- 
ous heterostructures [I?} , but these will not be discussed 
in this paper. Section 4 provides a summary and outlook. 
In the Appendix the discretization of the bath Green's 
function is discussed in more detail for the degenerate 
three-band model and the square lattice. Energy units 
are eV unless noted otherwise. 



ED DMFT FORMALISM 



2.1. Multi-band Correlations 



In this section we specify the details of ED DMFT 
appropriate for multi-band materials. In particular, we 
consider transition metal oxides, such as V2O3, LaTi04, 
and Ca 2 Ru04, where, as a result of crystal field split- 
ting, Coulomb correlations take place primarily within 
the partially occupied i 2 g bands. In addition to these 
three-band systems, we also consider the more general 
case in which all five d bands are relevant near the Fermi 
level. This applies to ordinary transition metals, oxides 
such as MnO or LaMn03, and iron pnictides and chalco- 
genides, such as FeAsLaO and FeSe. We focus on the in- 
fluence of onsite Coulomb correlations and neglect spatial 
fluctuations. Within this multi-band single-site DMFT 
approach the complex self-energy accounts for dynami- 
cal correlation effects, such as spectral weight transfer 
between low-energy and high-energy regions and quasi- 
particle broadening, while its momentum dependence is 
neglected. 

The single-particle properties are assumed to be given 
in terms of an effective tight-binding Hamiltonian i(k) 
which is derived from a fit to results obtained within 
density functional theory (DFT). The interacting Hamil- 
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tonian is given by: 

H = ^ t mn (]&)c£ mir Ck ncr -f- / ] Uf)'im\'ft'im\. 
mrik.tr im 

imer im<m' aa' 

(i) 

where c^ CT are annihilation (creation) operators for elec- 
trons on site « in orbital m with spin cr and ni mCT = 
c iLcr c inio" c kmcr are tne corresponding Fourier compo- 
nents, /i is the chemical potential and the intra-orbital 
and inter-orbital Coulomb energies are denoted by U 
and V, respectively. The exchange integral is J, where 
U' = U — 2 J because of rotational invariance [38| . 

For simplicity we limit the discussion to systems in 
which non-diagonal components of the density of states 
vanish for symmetry reasons. The diagonal elements of 
the lattice Green's function are then given by 

G m (iu] n ) = ^2[iu) n + M - t(k) - E(iw n )]m m , (2) 

k 

where w„ = (2n + l)ir/{3 are Matsubara frequencies 
(/3 = 1/T), and E is the self-energy matrix. Since local 
Coulomb interactions preserve the symmetry properties 
of G, E is also diagonal in orbital space, with elements 
E TO (iw n ). We consider only paramagnetic systems, so 
that the spin index of G and E can be omitted. Note that, 
because of the non-diagonal nature of t(k), the compo- 
nents of G are influenced by all components of E. In the 
limit of large ui n , the lattice Green's function components 
take the simple form 

G m {iLO n ) = [ioj n + n - e m - Y, m {iui n )]~ x , (3) 

where e m = E k t(k)] mm . 

For the purpose of the quantum impurity calculation 
within DMFT it is necessary to remove the self-energy 
from the central site in order to avoid double-counting 
of Coulomb interactions. This step yields the so-called 
bath Green's function 39] 

G , m (iu n ) = [G m (iw„) _1 + E m (iw n )] -1 . (4) 

Note that Go is not the non-interacting Green's function 
obtained from Eq. ([2]) for E = 0. Instead, Go accounts 
for the electronic motion from the central site throughout 
the lattice back to the central site, where E is present at 
all sites except the central one. 

Within ED DMFT the infinite lattice surrounding the 
central site is approximated via a bath of finite size. 
Thus, the above Green's function is approximated in 
terms of an Anderson impurity model for a cluster con- 
sisting of n c impurity levels £ m =i,...,n c and rc.& bath lev- 
els £fc=n c +i....,n s , which hybridize via hopping matrix el- 
ements V m k- Here, n s = n c + ni, defines the total number 




FIG. 1: (Color online) Left panel: Impurity bath interactions 
within diagonal basis for four orbitals. Green dots: impurity 
levels (m = 1, . . . , 4); blue dots: bath levels (k = 5, . . . , 12). 
Each impurity orbital interacts with its own bath via hopping 
terms V m k- The impurity orbitals interact only via Coulomb 
and exchange terms since the system is diagonal in this ba- 
sis. Right panel: Non-diagonal mixed basis involving impu- 
rity sites (i = 1, ... ,4) (red dots) and molecular-orbital bath 
levels (blue dots). Interactions among sites are determined 
by single-particle hopping, whereas the Coulomb interaction 
is purely onsite. Bath levels now hybridize with all impurity 
sites. The interactions indicated in the left and right panels 
are equivalent, since they differ only by a unitary transforma- 
tion among impurity orbitals. 



of cluster levels. The bath levels and hopping terms 
V m k may be viewed as auxiliary quantities to achieve an 
optimal representation of the lattice bath Green's func- 
tion in terms of a finite cluster. (The term "cluster" 
here refers to the finite set of impurity and bath levels 
and should not be confused with the "correlated impurity 
cluster" discussed in the next subsection.) Thus, 

G , m (iu> n ) « Gg' m (iw„), (5) 

where 

Go m (iu) n ) = (iui n + (i-e m - V* - — — — ) . (6) 

k—n c + l 

(The bath levels are defined relative to /u.) Since the 
lattice Green's function is diagonal in the orbital index, 
we can assume that each impurity orbital couples only 
to its own bath. For instance, for n c — 4 and rib = 8 
(i.e., two bath levels per impurity orbital), only the ele- 
ments Vm^m+A, Vm,m+8 7 an d their transpose elements are 
nonzero. All other elements vanish. These interactions 
are schematically indicated in the left panel of Figure 1. 
Details concerning the discretization indicated in Eq. ([6]) 
will be discussed in Section 2.4. 

Adding the onsite Coulomb interactions to the single- 
particle Hamiltonian specified by e OT , and V m ki the 
interacting cluster Green's function at finite T is derived 
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from the expression 2^, 21 1 

G d (luJ ) = L^ e -PB„( 1(m|c+»I 

/'■ 

\(n\c m aW)\ 2 



— E v 



E v — E,, 



(7) 



iGt+{^n) + Gt-{iu n )], 



where E v and \i>) denote the eigenvalues and eigenvec- 
tors of the full Hamiltonian, and Z = exp(— f}E v ) is 
the partition function. G^ v {iuj n ) are Green's functions 
derived from the excited states \v). G cl satisfies the same 
symmetry properties as G and Go, so that it is also di- 
agonal in orbital space. The orbital components of the 
cluster self-energy can be obtained from an expression 
analogous to Eq. (J4j) : 

E£(iw n ) = GS| m (iw n )- 1 - G^(kj n )-\ (8) 

The key assumption in DMFT is that this cluster self- 
energy is a physically reasonable representation of the 
lattice self-energy. Thus, 



E m (iw n ) = E^(iw n ), 
so that, as a result of Eq. ([3]), 

G m (iuj n ) » G^(iw n ). 



(9) 



(10) 



In the next iteration step, these diagonal self-energy com- 
ponents are used as input in the lattice Green's function, 
Eq. ©. The iteration cycle then has the schematic form 



Gfji-j Go,m 



Ill) 



In practice, the initial cluster parameters e m , and 
T^tj/b are obtained from a fit of the uncorrelated lattice 
Green's function or from a previous solution for weaker 
interaction parameters, other doping concentrations, etc. 

The ED DMFT self-consistency procedure involves at 
each iteration two projections: from lattice to cluster, 
and back from cluster to lattice. By definition, lattice 
quantities such as G, Go, and £ have continuous spectra 
at real energies, while the analogous cluster quantities, 
G cl , Gq, and E c/ are discrete. In this sense, the cluster 
quantities are never equal to the corresponding lattice 
versions, except for an infinite number of bath levels. 
The advantage of working at finite T is that, along the 
Matsubara axis at not too low temperatures, both sets 
of quantities become rather smooth, so that they can be 
meaningfully compared, even if only a finite number of 
bath levels is employed. 

Mathematically, the lattice-to-cluster and cluster-to- 
lattice projections in Eqs. §5§ and © imply that contin- 
uous and discrete real-energy spectra exist which have 



nearly identical representations along the Matsubara 
axis. In particular, the success of ED DMFT as a method 
for describing strongly correlated properties depends cru- 
cially on the accuracy of the discretization of Go , Eq. ([5]) , 
which, in turn, depends on bath size and temperature. 
This will be discussed in more detail in Section 2.4. The 
quality of this discretization will also be illustrated in 
Section 3 for a variety of realistic systems. Additional de- 
tails concerning the finite-temperature multi-orbital ED 
DMFT approach can be found in Ref. |2(| ■ 

Comparing the ED DMFT formalism outlined above 
to QMC DMFT, it is clear that both schemes involve 
additional errors beyond the basic single-site DMFT ap- 
proximation: In the case of ED, the bath Green's func- 
tion Go is discrctized, so that the interacting impurity 
Green's function suffers from finite-size effects. In the 
case of QMC, the discretization is avoided, but the in- 
teracting impurity Green's function suffers from errors 
due to statistical uncertainties. As a result, the limit 
of very low temperatures and energies is inaccessible to 
both schemes. 

We also note that the iterative procedure indicated in 
Eq. (|12p is not unique. For instance, in the definition of 
the cluster self-energy, Eq. ([8]), one could replace Gg m 
by the more accurate lattice function Go,m- Neverthe- 
less, Eq. ((HJ is found to yield a more reliable convergence 
towards a self-consistent solution. In addition, instead 
of finding the cluster parameters e m , and V m k m 
the (i + l) th iteration step from the actual lattice bath 
Green's function Go, m as indicated in Eq. ([5]), it is more 
convenient to employ an admixture 



G 



where Gq'^ is the cluster bath Green's function of the 
previous iteration step and a w 0.5. With these defini- 
tions, 5 to 10 iterations are usually sufficient to achieve 
convergence. 

As mentioned above, we have assumed so far that the 
local Green's function and self-energy arc diagonal in the 
orbital basis. Complex lattice geometries, however, can 
give rise to non-diagonal components which require a 
suitable generalization of the above formalism. This will 
be briefly addressed at the end of the next subsection. 



ci,0'+l) 
0,m 



aG , m + (1 - a)G 0> 



(12) 



2.2. Multi-Site Correlations 

Let us consider now a single-band Hubbard model 
and include spatial correlations among a finite number 
of sites (n c ). To be specific, we limit the discussion to 
two-dimensional systems, such as the square lattice with 
isotropic nearest and next-nearest neighbor coupling, and 
the triangular and honeycomb lattices with isotropic 
nearest-neighbor hopping integrals. The Coulomb inter- 
action within the band is assumed to be a purely on-site 
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manner, with elements E m . In the molecular-orbital ba- 
sis, the lattice Green's function takes the form: 



FIG. 2: (Color online) Unit cells within CDMFT for triangu- 
lar, square, and honeycomb lattices. 



interaction. Models of this kind are believed to capture 
important features related to short-range correlations in 
high-T c cuprate superconductors, organic molecular crys- 
tals, and in graphene. The Hamiltonian of the system is 
given by 

H = - Uj {cf^Cjcr + H.c.) + U ^ n it n n - X! V n icr 

(ij)a i i<r 

(13) 

where the sum in the hrst term extends up to first or 
second neighbors. 

Within the cluster extension ^ of DMFT (CDMFT) 
the interacting lattice Green's function in the cluster site 
basis may be written as 



G i:i (ioj n ) = J2i iuJ n + M - - E(iw„)]y\ (14) 

k 



where the k' sum covers the reduced Brillouin Zone and 
i(k') denotes the hopping matrix for the superlattice. £y 
represents the cluster self-energy matrix. The site labels 
for the three lattices are indicated in Figure [2] For the 
purpose of the ED calculation, it is convenient to change 
from the non-diagonal site basis to a "molecular-orbital" 
basis, in which the lattice Green's function is diagonal. 
For the triangular lattice, the diagonal components G m 
are given by: 

G\ = Gn + 2C?i2 
G2 = G3 = G11 — G13. (15) 

For the square lattice, the diagonal components are: 

Gi ; 2 = G11 i 2G12 + G14 
G3 = G4 = G11 — G14. (16) 

For the honeycomb lattice, they are given by: 

Gi,2 = (G11 + 2G13) ± (G14 + 2G12) 
G3,4 = G5.6 = (G11 — G13) ± (G14 — G12). (17) 

The cluster self-energy can be diagonalized in the same 



G m (iw„) = 2_J[ itjJ n + M - - E(iw„) 



-1 

mm 5 



(18) 



where i'(k') = T" 1 i(k')T and T = [T im ] denotes the uni- 
tary transformation from the site basis to the molecular- 
orbital basis. For the triangular, square and honeycomb 
lattices, the matrices T and the molecular-orbital com- 
ponents p m (w) of the uncorrelated density of states can 
be found in Refs. The asymptotic behavior of 

the Green's function components is given by 



G m (iuj n ) = [iu n + [i 



Srr.K)]" 1 , (19) 



where e m = E k < £'(k')] mm . 

Evidently, with the transformation from lattice sites 
to molecular orbitals, the ED CDMFT procedure be- 
comes analogous to the multi-orbital formalism discussed 
in the previous subsection. In particular, the removal 
of the self-energy from the "impurity cluster" yields the 
bath Green's function components Go, m («w rl ) given by 
Eq. (H]), which are discretized according to Eq. (JB]) by 
replacing the host continuum by a finite set of levels 
Thus, cluster molecular orbitals and bath form a "super- 
cluster" of size n s = n c + nb whose many-body properties 
can be solved exactly. (In the following we use the term 
"cluster" mainly for the unit cells depicted in Figure 2. 
From the context it will be clear whether the larger su- 
percluster is implied instead.) 

Assuming independent baths for the impurity molecu- 
lar orbitals, the hybridization interaction takes the form 
as schematically indicated in the left diagram of Figure 1 
for the case n c = 4. The advantage of this diagonal ba- 
sis is that, as in the multi-orbital case discussed in the 
previous subsection, the bath Green's function compo- 
nents GQ m (iuj n ) can be discretized independently, rather 
than fitting the non-diagonal Green's function matrix 
(iuj n ) in the site basis in one step (30L Eoj]. 
The evaluation of the interacting cluster Green's func- 
tion G cl can be carried out within the diagonal molecular- 
orbital basis, yielding the components G^iuin) defined 
in Eq. ([5]). In this basis, the onsite Coulomb interaction 
is a non-diagonal matrix, with intra-orbital and inter- 
orbital elements that are analogous to the Coulomb and 
exchange interactions of the single-site, multi-orbital sys- 
tems discussed in the previous subsection. For instance, 
in the case of a 2 x 2 cluster representing a square lattice, 
one finds (toiTO2||'ti3'714) = U/A for 64 of the 256 matrix 
elements, where U — (H\\ii) is the on-site Coulomb inter- 
action. All other elements vanish. An equivalent scheme 
consists in transforming the cluster molecular orbitals 
back to the site basis, by using the transformation T in- 
troduced above. Thus, the bath molecular-orbital levels 
remain unchanged and the Coulomb interaction is purely 
onsite. In this mixed site molecular-orbital picture the 
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non-diagonal cluster Green's function is evaluated from 
the expression 



1 -pe u / Hf^K/^V) 

„+ 



{ v \4a\v)^\ C j°\ v ) 
ILJn + E u - E n 



(20) 



where Gy v {iw n ) are Green's functions corresponding the 
excited state Since G?j satisfies the same symmetry 
properties as Gij, it is diagonal within the molecular- 
orbital basis, whose elements G^ coincide with those de- 
rived via Eq. (jSJ). Within the non-diagonal site basis, 
the cluster bath hybridization takes the form indicated 
in the right diagram of Figure 1. The molecular-orbital 
bath levels do not mix within this representation, but hy- 
bridize with all cluster sites. Further details concerning 
the finite-temperature multi-site ED DMFT approach are 
provided in Ref. [41 1. 

The lattice examples considered above are special in 
the sense that a diagonal molecular-orbital basis exists. 
More general cases, such as the non-isotropic triangu- 
lar or square lattices, require suitable generalizations of 
the above formalism. As long as direct coupling between 
bath levels is excluded, the most general hybridization 
between impurity and bath levels is still represented by 
the right-hand diagram in Figure 1, but the hopping 
terms Vik obey weaker symmetry relations than before. 
As a result, the back-transformation to the molecular- 
orbital basis does not yield a diagonal Green's function. 
Because of the off-diagonal elements, the diagonal com- 
ponents can no longer be fitted independently of one an- 
other. Instead, the full bath Green's function matrix 
(or a suitable subblock) must be discretized in a single 
sweep [13, Ei-tlll El S3j - Similarly, complex lattice ge- 
ometries in multi-band single-site cases can give rise to 
non-diagonal inter-orbital components. Since the impu- 
rity bath hybridization is then also of the form as indi- 
cated in the right diagram of Figure 1, the discretization 
must be carried out in one step, rather than using inde- 
pendent fits of the Go, m as in the diagonal case. 



2.3. Size of Hilbert Space 

An important criterion for the accuracy of ED DMFT 
is the level spacing related to the finite size of the clus- 
ter consisting of impurity and bath. The dimension 
of the Hilbert space for n s = n c + rib cluster levels is 
2 2 ™ s . Since different spin sectors of the Hamiltonian do 
not couple, the largest sector for n s = 9 has dimension 
N = (9!/(4!5!)) 2 = 15876. For n s = 12, the dimen- 
sion is (12!/(6!6!)) 2 = 853776, while for n s = 15, it is 



(15!/(7!8!)) 2 = 41409225. Fortunately, the cluster Hamil- 
tonian is extremely sparse, with typically no more than 
20 to 30 finite off-diagonal elements per row. Thus, using 
32 parallel processors for n s — 15, storage requirements 
per processor are about the same as for full matrix di- 
agonalization for n s = 8 levels (N = 4900). Moreover, 
because of the Boltzmann factor in the cluster Green's 
function, Eqs. (JSJ and ([2T]) . only a limited number of 
excited states is required at finite T. The quantum im- 
purity calculation can therefore be carried out very effi- 
ciently by using the Arnoldi algorithm 33j, which gen- 



erates a finite number of eigenvalues and eigenvectors of 
large sparse matrices. The contributions to the cluster 
Green's function are then given by a superposition of 
excited state Green's functions G^ u (ioj n ) or Gf^(iu n ) 
which can readily be calculated by applying the Lanczos 
method. 

The key point here is that for n s — 9, . . . , 15 the com- 
ponents of the interacting cluster Green's function, if 
evaluated at real ui, contain vastly more poles than their 
non- interacting counterparts, G$ m . Thus, even if a par- 
ticular impurity level e m couples only to two bath lev- 
els Eis, the indirect coupling to the remaining baths via 
internal Coulomb and exchange interactions within the 
impurity ensures a nearly continuous spectral distribu- 
tion. In the case of multi-site correlations, the diagonal 
molecular-orbital basis implies many intra-orbital and 
inter-orbital Coulomb and exchange-like terms, whereas 
in the non-diagonal site basis the coupling between baths 
stems from the single-particle hopping between cluster 
sites. 

Dynamical mean field theory may be viewed as a 
scheme for evaluating the complex self-energy, which de- 
fines the redistribution of spectral weight arising from 
Coulomb interactions. According to Eq. (JSJ, the cluster 
self-energy is governed by the spectral details of the inter- 
acting cluster Green's function. Thus, for n s = 9, . . . , 15 
the nearly continuous spectral distributions of the com- 
ponents G^ ensure accordingly smooth distributions of 
Moreover, because of the high quality of these func- 
tions along the Matsubara axis, analytic continuation 
to real oj close to the Fermi energy is usually found to 
be rather stable. In principle, the cluster Green's func- 
tion and self-energy can be evaluated close to the real- 
energy axis which is a crucial advantage compared to 
QMC DMFT. On the other hand, as pointed out above, 
at each iteration in ED DMFT the cluster self-energy 
(with its discrete level spectrum) is assumed to provide 
a physically reasonable representation of the continuous 
lattice self-energy, in the sense that along the Matsubara 
axis both are nearly identical. Thus, as in QMC DMFT, 
analytic continuation of the self-energy is desirable and 
can be carried out using various methods. (A detailed 
discussion of the analytic continuation of the self-energy 
is given in Ref. 



45|.) Since the single-particle properties 



are known at real energies, no further approximations 
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are then required. This should be contrasted to analytic 
continuation of the lattice Green's function, which im- 
plies continuation of both self-energy and single-particle 
features. (The limit of weak interactions makes the dif- 
ference between these two methods particularly striking.) 

For some of the examples discussed in the next section, 
we have employed the routine ratint [4(i | for the extrap- 
olation of the self-energy and lattice Green's function to 
real ui. Typically, up to a few hundred Matsubara points 
are used, with a small frequency-dependent broadening. 
This procedure usually yields reliable results at low en- 
ergies, while considerable uncertainties arise in the inco- 
herent region at larger energies. 

The multi-orbital and multi-site scenarios discussed 
above differ greatly from the simpler single-band, single- 
site case. Even if at some finite temperature the lattice 
bath Green's function can be discretized accurately via a 
non-interacting cluster Green's function using few bath 
levels, the small cluster size yields a very coarse level 
spectrum. It is therefore essential to include more bath 
levels so that the larger Hilbert space guaranties better 
energy resolution. Nevertheless, it has been shown that 
integrated quantities, such as the critical Coulomb ener- 
gies U c i,c2 determining the T/U phase diagram for the 
Mott transition at half-filling, are already well converged 
using only 3 or 4 bath levels [l9[ (see also 21|). 



2.4. Bath Discretization 

The second important criterion for the accuracy of ED 
DMFT concerns the projection of the bath Green's func- 
tion components Go,m(*w n ) onto to a finite cluster, as 
indicated by the discretization in Eqs. and ((BJ. Ev- 
idently, the quality of this discretization depends sensi- 
tively on temperature and bath size. Along the Matsub- 
ara axis the spectrum of Go tTn (u>) is effectively sampled 
with Lorentzian functions of width uj n — (2n + X)ttT: 



lmG . m (iuj n ) — — [ dui- 

7T J L 



■ImG , m (w). (21) 



Thus, choosing a very low temperature gives rise to fine 
structure in Go m (iw n ) that can be fitted only with a 
sufficiently large number of bath levels. For the descrip- 
tion of multi-band or multi-site correlations with no more 
than two bath levels per impurity orbital it is therefore 
important to employ a Matsubara grid for not too low T. 
As will be discussed in more detail in the next section, 
applications to a variety of multi-band and multi-site cor- 
relation problems suggest that a temperature in the range 
T « 5, . . . , 20 meV (for typical band widths of the order 
of W = 2 eV) permits discretizations of reasonable ac- 
curacy. (Note that at very small uj n a sufficiently large 
number of k points must be used to get accurate lattice 
Green's functions, Eqs. © and ([l4"|). See also Eq. flU).) 



To determine the bath levels Sk and hopping terms 
V m k we minimize the difference 



M 



0,; 



1W„ 



(22) 



n=0 



where M w 2 10 is the total number of Matsubara points 
and the weight function W n = l/cJ n is introduced to 
give more weight to the low- frequency region j2(| . Since 
G^2i l/i^n at large u> n , it is appropriate to optimize 
the low-energy region. To carry out the above fit, we use 
the routine minimize provided in Ref. To start the 
iterative procedure, we use bath parameters obtained for 
the uncorrelated system, or from a converged solution for 
nearby Coulomb and exchange energies. The resulting Sk 
and V m k are usually very stable against (not too large) 
variations of initial conditions. 

We have also tested other weight functions, such as 
W n — 1 and W n = 1 /u) 2 , as well as minimization of [26[ 

M 

Dff4 = Y, WnlGoM^n)- 1 - G^ioJn)- 1 ] 2 , (23) 

n=0 

which represents the discretization of the lattice hy- 
bridization function 



A m (iw n ) = iu n + fjt - e m - G 0}m (ico n y 



(24) 



via the corresponding cluster expression, i.e., the last 
term in Eq. ©. Thus, 



A m (iuj n ) « A m (iuj n ) = V ,' mfc| , (25) 



so that 



M 

E 

71=0 



DifC = V W n \A m {tu n ) - A m {tu n )\ 2 . (26) 



It can easily be shown that, in the limit of large w„, 

\G 0>m (iu n ) - Gt m (iuj n )\ -t — \A m (iuj n ) - A m (iu n )\. 

(27) 

Thus, minimization of \G — G cl \ ensures greater accuracy 
at small io n than minimization of \1/G — 1/G |. Ex- 
perience has shown that the ED DMFT iteration proce- 
dure, using the distance function defined in Eq. (|22j) with 
W n = l/w„, yields reliable and efficient convergence. 
Nonetheless, as discussed in the Appendix, other weight 
functions, as well as utilization of Eq. (1231) . also pro- 
vide accurate discretization of the bath Green's function. 
(For a recent discussion of various choices of weight func- 
tions, mainly for the one-dimensional Hubbard model, 
see Ref. 

In principle, the asymptotic behavior of GQ. m (ioj n ) 
and Gq m (iujn) is governed by the impurity levels e m , 
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which follow from the single-particle properties of the 
multi-band or multi-site Hamiltonians, as indicated in 
Eqs. ^ and (fH?1) . Since this asymptotic behavior is only 
weakly affected by the precise value of e m , a more ac- 
curate discretization of Go,m(*w ra ) in the low-frequency 
region may be achieved by allowing e m to vary. In ef- 
fect, e m is then included among the auxiliary parame- 
ters used to obtain the optimal fit of the lattice bath 
Green's function Go, m (iu> n ) via the cluster bath Green's 
function Gq m (iui n ), even though this implies a (small) 
inconsistency between cluster and lattice properties. As 
a consequence of the finite bath size, slight inconsisten- 
cies of this kind are an intrinsic feature of ED quantum 
impurity calculations. They appear, for example, in the 
orbital occupancies, as calculated from the cluster and 
lattice Green's functions. In most cases we have found 
the freely determined impurity levels to be very close to 
the nominal bulk values specified in Eqs. ([3]) and (|T9| . 
Thus, the variation of e m is not essential for the dis- 
cretization procedure. 

For a typical situation with two bath levels per im- 
purity orbital, each complex function Go im (iui n ) is then 
fitted with five parameters: two bath levels two hop- 
ping terms V m k, and e m . For instance, in the case of 2 x 2 
clusters representing the square lattice with n c = 4 and 
nt, = 8, we have a total of 15 fit parameters for the three 
independent Green's function components (12 fit param- 
eters for fixed e rn ). Naturally, this parameter set provides 
a more accurate bath discretization than in the site ba- 
sis, where for symmetry reasons only four fit parameters 
are usually taken into account. In the special case of 
particle-hole symmetry at half-filling, the diagonal basis 
employs 8 to 10 parameters, compared to 2 parameters in 
the site basis, unless additional hybridization terms are 
incorporated as indicated in the right panel of Figure 1. 

3. RESULTS AND DISCUSSION 
3.1. Degenerate Three-Band Model 

As a first application of multi-band ED DMFT we con- 
sider a three-band Hubbard model consisting of identical 
subbands with a Bethe lattice density of states (t = 1/2, 
band width W — 2). The subbands are coupled only 
via Coulomb and exchange interactions. Single-particle 
hybridization among bands is neglected. Despite its sim- 
plicity, this model is highly non-trivial. As shown by 
Werner et al. [47], it exhibits metal insulator transi- 
tions at integer occupancies and a spin freezing transition 
within an intermediate range of occupancies. At the lat- 
ter transition, the self-energy changes from Fermi-liquid 
to non-Fermi-liquid behavior and the spin-spin correla- 
tion function changes from Pauli to Curie- Weiss behav- 
ior. We discuss this model here in some detail in order 
to demonstrate the convergence of the ED results with 





1 1 
ED 


I 1 

U = 4 T = 


0.01 






n = 1 .54 








.... 1.8 








2.0 


2.37 








2.6 



0.5 1 1.5 2 



\ QMC 


U = 


4 T = 0.01 


V 


n = 


1.6 


v ■**-.. 




1.79 






1.99 






--'--':::1M5—-. 


i 


i 


i 


3 0.5 


1 


1.5 









FIG. 3: (Color online) Self-energy of degenerate three-band 
model for U = 4, J = (7/6, T = 0.01 at several occupancies. 
Upper panel: ED DMFT results for n s = 12 (3 bath levels 
per impurity orbital). Lower Panel: continuous-time QMC 
DMFT results for similar occupancies (adapted from Figure 3 
of 0). 

bath size, the quality of the discretization of the bath 
Green's function Go, and the consistency with analogous 
continuous-time QMC results reported in Ref. [13] • We 
also discuss briefly that, for sufficiently large Hund's rule 
coupling, the self-energy reveals a collective mode due to 
spin fluctuations, which leads to a pseudogap in the den- 
sity of states. A similar collective mode is found in the 
five-band model discussed in Section 3.3. 

Figure [3] (upper panel) shows the subband self-energy 
as a function of Matsubara frequency for various occu- 
pancies. The cluster size is n s = 12, with three bath 
levels per orbital. Full Hund exchange is used, with 
J = £7/6, and the Coulomb energy £7 = 2W — 4 is 
chosen so that at half-filling the system is a Mott insu- 
lator (£7 C w 1.5 [48|). The corresponding QMC results 
are given in the lower panel for similar occupancies. The 
comparison shows that there is nearly quantitative agree- 
ment between the ED and QMC results, in particular, in 
the important low-energy region. At n rs 1.6, the system 
is a Fermi liquid, with ImS m (iw„) — > oj n in the limit of 
small u> n . Near n « 2, ImS,„ begins to exhibit a finite 
onset as a result of local- moment formation. At larger oc- 
cupancies, this onset continues to increase until the Mott 
phase is reached and, at half-filling, ImS m (zo; n ) — > l/ui n 
at small w„. The results in Figure [3] are for T=0.01. 
Similarly good agreement is obtained for T = 0.005 (not 
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FIG. 4: (Color online) Convergence of self-energy for three- 
band model with bath size at different chemical potentials. 
Solid red curves: n s — 9, dashed blue curves: n 3 = 12, long- 
dashed green curves: n a — 15, corresponding to two, three, 
and four bath levels per impurity orbital, respectively. For 
clarity, the results in the upper panel for ImS m at fi — 4 are 
shifted down by 0.5. U = 4, J = (7/6, T = 0.01. 



FIG. 5: (Color online) Discretization of bath Green's function: 
lattice Go,m(*Wn) (solid red curves) and cluster Go, m (* w n) 
(dashed blue curves) for 2 bath levels per impurity orbital 
(n s = 9) at several chemical potentials; U — 4, J = (7/6, 
T — 0.01. The dashed (magenta and green) curves show the 
fit for y, = 2 at T = 0.0025. (For clarity the latter results are 
shifted vertically by ±0.05). 



shown here). 

Analogous results for only two bath levels per orbital 
(n s = 9) show qualitatively the same behavior as for 
n s = 12. In fact, at fi = 2 (n ~ 1.6) and fi = 4 
(n w 2.5), the self-energies are nearly the same as for 
the larger cluster. Only in the transition region near 
fi = 3 (n « 2.0), larger differences are found, including 
a stronger dependence on temperature. To illustrate the 
convergence of the self-energy with bath size more sys- 
tematically, we show in Figure [4] results of ED DMFT 
calculations for n s = 9, 12, and 15. The calculations for 
2 and 3 bath levels per orbital are carried out at T = 0.01, 
whereas for computational reasons the ones for n s = 15 
are performed at T = in the evaluation of the cluster 
Green's function, Eq. <j8j> . The Matsubara grid employed 
in the DMFT self-consistency procedure corresponds to 
T = 0.01 in all cases. For = 2 and 4, both the real 
and imaginary parts of E m are in good agreement for 
the different cluster sizes. In particular, in the Fermi liq- 
uid range for fi = 2, ImS m ~ ui n in the limit of small 
U! n , with an effective mass enhancement of m*/m ~ 3. 
In contrast, for (j, = 4, ImE m approaches a constant, re- 
flecting the finite scattering rate due to the presence of 
frozen moments. For fj, = 3, the increase from n s = 9 to 
15 is seen to shift the minimum of ImE m (iu;„) to larger 



u) n and to make it more shallow, with a larger change 
from n s = 9 to 12 than from n s — 12 to 15. The reason 
for the stronger variation with bath size at fi = 3 is that 
this chemical potential nearly coincides with the phase 
boundary of the spin freezing transition. Thus, small 
variations due to the finite cluster size can give rise to 
larger changes of the self-energy. 

We point out that the self-energy shown in Figure 0] 
as a function of bath size does not exhibit any oscilla- 
tory behavior related to the even or odd number of bath 
levels per impurity orbital. This kind of dependency on 
bath size occurs in more special cases, such as particle- 
hole symmetric bands at half-filling. (See, for example, 
the two-band model discussed in Ref. |19|.) The pres- 
ence or absence of bath levels at the Fermi energy then 
tends to favor metallic or insulating behavior, leading to 
a slow convergence with the number of bath levels. In 
the present degenerate three-band model, the bath levels 
away from half-filling do not obey any particular symme- 
try rules, so that oscillatory variation with bath size does 
not arise. 

As discussed in the previous section, a crucial crite- 
rion for the convergence with bath size in ED DMFT is 
the quality of the discretization of the bath Green's func- 
tion Go tm (iu! n ), Eq. ([3]). via the cluster Green's function 
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FIG. 6: (Color online) Low-energy region of self-energy of 
three-band model for several chemical potentials at tempera- 
tures T = 0.01 (red solid curves), T = 0.005 (short-dashed 
blue curves), and T — 0.0025 (long-dashed green curves); 
U = 4, J — U/6, n a —9 (two bath levels per impurity or- 
bital). For clarity, the curves in the upper panel for y, = 4 are 
shifted down by 0.05. 



GQ m (iu> n ), Eq. (JSJ. These functions are compared in Fig- 
ure [5] for several chemical potentials, assuming two bath 
levels per orbital. Evidently, at T = 0.01, the fit of the 
bath Green's function with 5 parameters is already re- 
markably good. Including more bath levels improves the 
fit at low uj n even further. Also shown is the discretiza- 
tion of Go, m at T = 0.0025 for /i = 2. In this case the fit 
is much less satisfactory, indicating that more bath levels 
are required at this low temperature. 

The results in Figure 5 are based on minimization of 
\Gq — Gq\ as in Eq. (|2"2"|) . using the weight function W n = 
l/u>n \ where N = 1. In Appendix A.l it is shown that 
discretizations of similarly good quality are obtained for 
N = and N = 2, as well as for minimization of |l/Grj — 
1/G§'|, as in Eq. with N = 0,1, 2. 

It should be noted that the small kinks in Im Gg m (iu) n ) 
at low oj n are related to the fact that, because of the finite 
size of the cluster, the cluster spectrum is always gapped. 
Thus, ImGo TO (io;) — in the limit u — > 0, giving rise to 
large deviations from the lattice bath Green's function 
(see also Section 3.5 and Figure [21"]). These deviations 
occur primarily at very low frequencies near or below 
cjo- Naturally, they are much less pronounced if the sys- 
tem is insulating. To achieve an accurate discretization 
of Gq m with a small number of bath levels it is there- 




FIG. 7: (Color online) Spin-spin correlation function for de- 
generate three-band model at various occupancies (increasing 
from below) for n B = 12; U = 4, J = U/6, T = 0.01. In the 
Fermi-liquid regime at small n, the susceptibility is Pauli-like, 
whereas at larger n it reflects Curie- Weiss behavior. 



fore necessary to use a Matsubara grid for not too low 
temperatures. On the basis of the example shown in Fig- 
ure [S] and various other systems involving multi-orbital 
and multi-site correlations, T » 0.005, . . . , 0.01 is a rea- 
sonable lower limit, if two or three bath levels per orbital 
are available. 

Although the discretization of the bath Green's func- 
tion becomes less accurate at low temperatures, the self- 
consistently determined self-energy can still be qualita- 
tively correct. To illustrate this point, we show in Fig- 
ureEDthe self-energy at T = 0.01, 0.005, and 0.0025. Ev- 
idently, all curves exhibit the same trend as a function of 
chemical potential, with Fermi-liquid behavior at low oc- 
cupancy (/i = 2, n » 1.6) and a pronounced low-energy 
scattering rate closer to half-filling (// = 4, n « 2.5). 
As in Figure HI the results near the spin freezing phase 
boundary (fx = 3, n « 2) are most sensitive and therefore 
show larger variation with T. 

As shown in Ref. [47| . the physical reason for the 
change of the self-energy from Fermi-liquid to non- Fermi- 
liquid behavior seen in Figure |3] is associated with a 
qualitative change of the spin-spin correlation function 
G m (r) = (S mz (T)S mz (0)), where r denotes imaginary 
time. In the Fermi-liquid regime \x < 3, this function 
decays with r, so that G m (r = /3/2) is very small. The 
susceptibility Xm ~ J dr (5 m2 (r)S' mz (0)) then is Pauli- 
like, i.e., independent of temperature. For /i > 3, how- 
ever, G m (r) approaches a finite value for r — > /3/2, so 
that Xm becomes proportional to 1/T, as expected for 
Curie- Weiss behavior due to the formation of local mo- 
ments. Figure [7] shows the correlation function as cal- 
culated within ED DMFT. The variation of G m (r) with 
doping found here agrees rather well with the continuous- 
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time QMC results discussed in Ref. 

We note here that a similar sequence of phases, from 
a Mott insulator at half-filling towards non-Fermi-liquid 
and Fermi-liquid behavior with increasing electron or 
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Coulomb interactions in several materials with partially 
filled tig bands, where the e g subbands are pushed above 
the Fermi level as a result of crystal field effects. In this 
subsection, we summarize the main results of these calcu- 
lations, illustrate the convergence of the self-energy with 
the size of the bath, and point out the consistency with 
complementary QMC DMFT results. 

3.2.1. Ca 2 _ a; Sr I Ru04 



FIG. 8: (Color online) Imaginary part of self-energy of three- 
band model at U = 4, J = U/3 for several chemical potentials. 
T = 0.01, n s — 12. Note the onset and low-energy kink. 



hole doping, is also seen in the five-band Hubbard model 
discussed in Section 3.3 and in the non-local properties 
of the single-band Hubbard model discussed in Section 
3.5. 

Apart from the spin freezing transition, the three- 
band model exhibits another interesting phenomenon not 
present in the single-band case. As shown in Figure[8j for 
Hund exchange J = U/3 at U = 4, the imaginary part of 
the self-energy reveals not only a finite low-energy scat- 
tering rate, but also a kink at about cj n pa 0.2. Extrap- 
olation of the self-energy to real lu yields a resonance in 
ImE(u) near w = —0.2 (below Ep — 0), which is weak 
in the Fermi-liquid region (n < 2), but gets very pro- 
nounced in the non- Fermi- liquid range (n > 2). More- 
over, the resonance disappears at small J. The reason for 
this dependency on J is the fact that, at half-filling, the 
three-band Hubbard model is metallic at small J, and 
becomes Mott insulating for sufficiently large J 
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Thus the resonance exists only in proximity to the Mott 
phase in the limit of half-filling. It can also be shown 
that the peak in ImS(w) gives rise to a pseudogap in 
the density of states slightly below Ep and that it is as- 
sociated with a maximum in the local spin correlation 
function. The resonance in the self-energy induced via 
Hund coupling may therefore be viewed as a collective 
mode caused by spin fluctuations. According to Figure 
IHl the non- Fermi-liquid behavior associated with the low- 
energy onset of Im£(w) can then be interpreted as a 
consequence of the appearance of this collective mode. 
A similar phenomenon occurs in the five-band Hubbard 
model discussed in Section 3.3. For a discussion of col- 
lective modes due to spin fluctuations in the single-band 
Hubbard model, see Ref. (50l . 



3.2. Three-Band Materials 

The multi-band ED DMFT approach outlined in Sec- 
tion 2.1 has recently been used to study the role of local 



A system that has been studied extensively is the layer 
perovskite Ca2- 2; Sr 2 ;Ru04 which exhibits a remarkably 
rich phase diagram, ranging from unconventional super- 
conductivity in the pure Sr limit to an antiferromagnetic 
Mott insulator at x — 0, with a variety of exotic magnetic 
phases at intermediate concentrations [5l|. As the sub- 
stitution of Sr by the slightly smaller Ca ions is isoelec- 
tronic at occupancy 4d 4 , the various electronic and mag- 
netic phases are entirely the result of subtle hybridization 
changes among the ti g states caused by structural distor- 
tions. Because of the essentially two-dimensional struc- 
ture of this material, the d xy band has a much larger 
band width than the more nearly one-dimensional d xz>yz 
bands. Thus, Coulomb interactions can be expected to 
modify these subbands differently. Moreover, LDA cal- 
culations [H^] predict the d xy van Hove singularity for 
x = 2 to lie only about 50 meV above the Fermi level, 
so that correlation induced charge transfer between the 
tig bands could possibly shift this singularity below Ep, 
giving rise to a qualitative change of the topology of the 
d xy Fermi surface from electron-like to hole-like. Multi- 
band QMC and ED DMFT calculations for Sr 2 Ru0 4 do 
indeed predict a decrease of orbital pola rization, i.e., a 
shift of charge from d xz . yz to d xy [53145 71 ]. The concomi- 
tant lowering of the d xy van Hove singularity for realistic 
interaction parameters, however, is not large enough to 
modify the topology of the Fermi surface. 

The upper panel of Figure [9] illustrates the conver- 
gence of the ED self-energy components with cluster size 
for a tight-binding Hamiltonian including up to second- 
neighbor hopping interactions 53(. The results are for 



full Hund exchange and the values of U and J corre- 
spond to those derived in recent constrained RPA calcu- 
lations [I3|- There is good agreement for two and three 
bath levels per impurity orbital, especially at low ener- 
gies. The small differences are mainly caused by slightly 
different orbital occupancies obtained for these two bath 



sizes. For n s = 12, we find n n 



0.66 and n xz 



The uncorrelated subband occupancies are n x 



0.67. 

: 0.63 



and 



L xz,yz 



0.68. Thus, Coulomb correlations in the 



pure Sr compound lead to decreasing orbital polariza- 
tion. The effective mass enhancement is approximately 
m*/mLDA w^in agreement with continuous-time QMC 



calculations [57j and with experiment [51]. In this sense, 
it is appropriate to regard Sr2Ru04 as a strongly corre- 
lated material. 
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FIG. 9: (Color online) Upper panel: Imaginary part of 
Sr2Ru04 t2 g self-energy components for cluster sizes n s = 9 
and n s = 12; U = 2.3 eV, J = 0.4 eV, T = 5 meV. Lower 
panel: Orbital occupancies (per spin) of Ca2- a; Sr a; Ru04 ti g 
subbands as functions of Coulomb energy for various crystal 
field splittings A, derived within ED DMFT (n s = 9) [He|. 
Solid (red) dots: n xy , empty (blue) dots: n xz , y z. Hund ex- 
change J = [7/4; T = 20 meV. 



For a tight-binding model similar to the one used in 
Ref. [53j], continuous-time QMC results [HH for the self- 
energy agree very well with those in Figure In particu- 
lar, they show little difference between the d xy and d xz%yz 
effective mass enhancements. Stronger ti g anisotropy is 
obtained for a more accurate single-particle Hamiltonian, 
for which the d xy van Hove singularity lies closer to Ep 
and the difference in d xy and d xz . yz band widths is less 
pronounced than in the tight-binding model [571 ]. 

Substitution of Sr via Ca leads to rotation, tilting and 
flattening of Oxygen octahedra and, as a result, to a low- 



crm 



of the d xy band with respect to the d x 



bands 



59l 60j . To investigate the influence of Coulomb correla- 
tions in the presence of this structurally induced interor- 
bital charge transfer we have performed ED DMFT calcu- 
lations, with an additional crystal field A = e 



xz.yz 



- .ry 



to 



account for the downward shift of the d xy bands [56]. As 
shown in the lower panel of Figure |H1 correlations greatly 
enhance the d xx ,yz ~~ ^ d xy charge transfer (except at small 
values of U). In fact, for realistic values of A, U, and 
J, the d xy band gets completely filled and the half-filled 
d XZ:VZ bands become Mott insulating. This orbital po- 
larization is also evident in the cluster spectra which are 
depicted in FigurcHUlfor A = 0.4 eV. At U = 3.0 eV, the 



FIG. 10: (Color online) Quasiparticle spectra of Ca2Ru04 for 
crystal-field splitting A = 0.4 eV at U = 3.0 eV (upper panel) 
and U = 4.5 eV (lower panel) for J = (7/4, T = 20 meV. Solid 
(red) curves: d xy band, dashed (blue) curves: d X z,yz bands. 
The spectra are derived from the interacting cluster Green's 
function for n s — 9, with broadening 7 = 50 meV [H^ ]. 



system is metallic, and all three ti g orbitals contribute 
to the density of states at Ep. At U — 4.5 eV, how- 
ever, the d xy band is pushed below Ep and the doubly 
degenerate d xz . yz bands are split into lower and upper 
Hubbard bands. Within the accuracy of the ED DMFT 
calculations, the filling and Mott-splitting of subbands 
takes place at the same critical U c ~ 4.2 eV. Thus, there 
are no orbital-selective, successive Mott transitions for 
the narrow d xz . yz and wide d xy subbands, in contrast to 
the mechanism proposed by Anisimov et al. ,54]. Re- 
cent continuous-time QMC DMFT results for Ca2Ru04 
by Gorelov et al. [13] based on detailed LDA calculations 
also do not find successive transitions. Orbital selective 
Mott phases can appear, however, if the t 2g crystal field 
splitting is assumed to have the opposite sign 61]. 

We point out that the large critical Coulomb energies 
seen in the lower panel of Figure [9] for small values of A 
are a consequence of the assumption J — U/4. If J is held 
fixed at a realistic value, U c becomes accordingly smaller. 
For U = 3.1 eV and J = 0.7 eV, as in Ref. [6<|, we have 
J r; [7/4 so that the orbital polarization obtained for 
A w 0.4, . . . , 0.5 eV is qualitatively correct. 



3.2.2. LaTi0 3 

The nature of the Mott transition in Ca2Ru04, where 
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FIG. 11: (Color online) Upper panel: LaTiOa subband oc- 
cupancies for orthorhombic and tretragonal bulk geometry. 
Since the crystal field splitting has the opposite sign in these 
structures, the orbital polarization between the singly degen- 
erate a g (dx V ) and doubly degenerate e g i (d X z,yz) states is 
reversed [63] ■ Lower panel: LaTiOa *2 9 subband occupancies 
per spin as functions of chemical potential. (7 = 5 eV, Hund 
exchange J = 0.65 eV, T — 20 meV. Also shown is the av- 
erage occupancy n av = (n ag + 2n e9 <)/3. The arrows indicate 
the hystereses for increasing and decreasing /x [6(| . 



the correlation induced crystal field enhancement leads to 
strong orbital polarization, is similar to the one discussed 
by Pavarini et al. 62J for LaTiC>3 (3d 1 ). In contrast to 
cubic SrVC>3, the orthorhombic structure of LaTiOa gives 
rise to a splitting of about 0.2 eV between the singly de- 
generate a g and doubly degenerate e g ' states, implying 
appreciably different subband occupancies: n ag sa 0.23 
and n e , ps 0.135 per spin. This orbital polarization 
is strongly enhanced by the local Coulomb interaction, 
leading to virtually empty e g > subbands and a Mott phase 
in the half-filled a g band. Finite-temperature ED DMFT 
results [63J which are shown in upper panel of Figure [TT] 
are consistent with this picture. 

Interestingly, when thin layers of LaTiOa are deposited 
on cubic SrTiC-3, as in LaTiOa/SrTiOa heterostructures, 



the initial growth appears to be tetragonal [6J|. This 
substrate-enforced reconstruction gives rise to an in- 
creased ti g band width and to a reversal of the ti g crys- 
tal field. As also shown in the upper panel of Fig- 
ure [Til the result of these effects is an increase of the 
critical Coulomb energy for the Mott transition, and 
qualitatively different subband occupancies. Now the 



singly-degenerate d xy band is emptied, and the quarter- 
filled d XZtVZ bands are split into lower and upper Hub- 
bard bands. These results are relevant for the interpre- 
tation of the metallicity observed LaTiOa/SrTiOs het- 



erostructures [65j . If realistic values of U for LaTiC-3 are 
about 5.0, . . . , 5.5 eV, the results in Figure [TT1 imply that 
thin LaTiOa layers are strongly correlated metals rather 
than Mott insulators. Thus, the experimentally observed 
metallic behavior is not associated with the atomic inter- 
face Ti layers with nominal 3d 5 occupancy, but with the 
entire LaTiOa slabs. 

The lower panel of Figure [TTJ demonstrates that the 
ti g orbital polarization is also highly sensitive to electron 
and hole doping [66| . Since the subband compressibili- 
ties K m = dn m /d[i diverge at the Mott transition, small 
changes away from the nominal n — 1 occupancy cause 
significantly larger changes in subband occupancies. In 
fact, weak overall hole doping gives rise to much larger 
hole doping in the half-filled a g band, accompanied by 
electron doping in the e g i bands. Conversely, weak over- 
all electron doping of the system gives rise to much larger 
electron doping in the empty e g > bands, accompanied by 
hole doping of the a g band. Thus, near the Mott tran- 
sition, pronounced internal charge rearrangement takes 
place among the relevant subbands. 

ED DMFT cluster spectra illustrating the doping de- 
pendence in LaTiOs, and cluster spectra corresponding 
to the orthorhombic and tetragonal bulk phases, can be 
found in Refs. [6r3 | and [63j |. respectively. 

3.2.3. V 2 3 

Another material that exhibits enhanced orbital po- 



larization near the Mott transition is V 2 3 (3d 2 ) [67|. 
In this case, the corundum lattice geometry ensures that 
the doubly-degenerate e g > bands have a slightly larger 
binding energy than the singly-degenerate a g band. In 
Ref. [lo], we have used ED DMFT to analyze this or- 
bital polarization for n s = 9, i.e., using two bath levels 
per impurity orbital. Figure [T2l illustrates the quality of 
the discretization of the bath Green's function Go tm (ico n ) 
components. As pointed out above, cluster spectra are 
always gapped, so that ImGg' m (zo; n ) — )• at very low 
frequencies. This limiting behavior is discernible, how- 
ever, only below ojq = irT = 0.031 eV, where the lattice 
and cluster Green's functions strongly deviate from one 
another. Thus, by choosing a Matsubara grid with a 
temperature not much lower than T w 10 meV, these 
discrepancies are almost completely avoided and the fi- 
nite cluster can serve as an adequate representation of 
the lattice continuum. 

Figure [13] shows the Coulomb driven enhancement of 
the V2O3 orbital polarization. Near U ~ 5 eV, the 
a g band is pushed above the Fermi level and the half- 
filled e g i bands are split into lower and upper Hubbard 
peaks. Note that the precise value of the critical Coulomb 
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FIG. 12: (Color online) Discretization of V2O3 bath Green's 
function components, using two bath levels per impurity or- 
bital; U = 5 eV, J = 0.7, T = 10 meV. Upper (lower) panel: 
imaginary (real) part. 



energy depends sensitively on the exchange terms in- 
cluded among the many-body interactions. In the case of 
only density-density terms (Ising-like exchange), we find 
U c ~ 5.1 eV, in agreement with Hirsch-Fye QMC DMFT 
calculations by Keller et al. 67]. Full Hund exchange, 
on the other hand, shifts U c upward to about 5.6 eV, a 
trend that is consistent with results obtained in ED and 
NRG studies of simpler two-band models 19l. 68|. 

The lower panel of Figure Q2] shows the convergence of 
the self-energy components in the metallic phase close to 
the Mott transition with bath size. The nearly half-filled 



e g i subbands are seen to exhibit a significant low-energy 
scattering rate which is reminiscent of the spin freezing 
transition in the degenerate three-band model discussed 
in Section 3.1. The nearly empty a g band, on the other 
hand, is much more itinerant. Nevertheless, it exhibits a 
large effective mass enhancement of about 5, and a small 
but finite scattering rate caused by the interaction with 
the almost Mott-localized e g > states. The differences for 
n s = 9 and n s — 12 are presumably caused by slightly 
different orbital occupancies. Close to a transition, the 
self-energy is particularly sensitive to small parameter 
changes. A similar effect was found close to the spin- 
freezing transition in the degenerate three-band model 
discussed in Section 3.1. 



3.2.4. Nao. 3 Co0 2 



FIG. 13: (Color online) Upper panel: V2O3 subband occu- 
pancies per spin as functions of Coulomb energy. Red sym- 
bols: Ising exchange, blue symbols: full Hund exch ang e. Solid 
dots: ED DMFT gcj, empty: dots: QMC DMFT [gTfl. Lower 
panel: imaginary part of self-energy components for two and 
three bath levels per impurity orbital (n s = 9 and 712 = 12, 
respectively) at U = 5 eV, J = 0.7, T = 10 meV. 



Despite several years of experimental and theoretical 
investigations, basic features of the electronic properties 
of the intercalated layer compound Na^OoO^ are still 
controversial. The doping concentration x = 0.3 is of 
particular interest since this material becomes supercon- 
ducting when small amounts of water are added [691 ] . 

As a result of the octahedral crystal field, the Co 3d 
bands are split into t 2g and e g subbands. With Na dop- 
ing the occupancy of the t2 g sector can be continuously 
tuned between n — 5 and n — 6. The rhombohedral 
symmetry at Co sites yields a further splitting of ti g or- 
bitals into singly-degenerate a g and doubly-degenerate 
e g > states. Because of the nearly two-dimensional hexag- 
onal structure, the Fermi surface as predicted by LDA 
calculations [70| exhibits a large a g hole pocket centered 
at r, and six small e g < hole pockets along the TK direc- 
tions of the Brillouin Zone. These small pockets have so 

0, Ill- 



far not been observed in ARPES measurements 
As discussed above, local Coulomb interactions can give 
rise to charge transfer between subbands, implying in- 
creasing or decreasing orbital polarization, depending on 
the details of the single-particle Hamiltonian. Several 
groups have carried out DMFT calculations in order to 
explore whether the small pockets might disappear be- 
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FIG. 14: (Color online) Upper panel: Nao.3Co02 subband 
occupancies as functions of Coulomb energy for two single- 
particle Hamiltonians Hi and H2 (see text) and different 
Hund exchange energies J (74|. Lower two panels: Real and 
imaginary components of Nao.3Co02 self-energy calculated 
within ED DMFT for two and three bath levels per impu- 
rity t2 g orbital; Hamiltonian H2, U = 3 eV, J — 0, and 
T = 20 meV. 



cause of a q to e q ' electron transfer 



The upper panel of Figure [T4l summarizes the degree of 
orbital polarization obtained for slightly different single- 
particle Hamiltonians Hi and H2, as well as different 
Hund couplings J 74]. The important difference be- 
tween these Hamiltonians is that, because of different 
band crossings in the tight-binding fits, the effective crys- 
tal field A = e a — e e . is about —130 meV for Hi, while it 

tig C g l ±1 

is —10 meV for Hi- The band crossings obtained for Hi 
correspond to those of high-quality LAPW LDA calcu- 
lations. As shown in the figure, these two tight-binding 
fits yield opposite trends of the orbital occupancy with 
increasing local Coulomb energy. Hi gives rise to de- 
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FIG. 15: (Color online) Quasiparticle spectra for a g (up- 
per panel) and e' g (lower panel) subbands of Nao.3Co02 at 
(7 = 3 eV, J = 0.75 eV, T = 10 meV. Solid (red) curves: spec- 
tra derived from lattice Green's function G m (w) after extrap- 
olation of E m («*;„) to real cj; dashed (blue) curves: spectra 
derived via extrapolation of G m (iu} n ); dotted (black) curves: 
single-particle density of states [2(| . 



creasing orbital polarization, whereas Hi leads to increas- 
ing orbital polarization. These different trends are main- 
tained even if the Coulomb and energies are varied over 
wide ranges, U m 3, . . . , 6 eV and J = 0, . . . , 0.9 eV. Ev- 
idently, correlation effects are not strong enough to fill 
the e g / pockets. 

The lower panels of Figure [TJ] illustrate the conver- 
gence of the ED self-energy components of NazCoC^ with 
cluster size, using Hamiltonian H^. The small differences 
obtained for n s — 9 and n s = 12 are partly related to the 
slightly different occupancies: n ag = 0.753 (0.740) and 
n dg = 0.946 (0.959) for n s = 9 {n s = 12). The impor- 
tant low-frequency features are seen to be in very good 
agreement for these cluster sizes, in particular, the small 
upward shift of a g states relative to e g i by about 0.15 eV. 
We also point out that these ED results for H 2 are fully 
consistent with the ones obtained within QMC for the 
same Hamiltonian j76[. In fact, the detailed comparison 
of the respective self-energy components shows, that, if 
identical input parameters are employed, nearly quanti- 
tative agreement is found for these impurity solvers 



77] 



adjusting the 
. As a result, 

the e g > bands are pushed down, and Coulomb correlations 
can fill them completely. The physical origin of such a 



The e g r pockets can be made to vanish by 
h g crystal field to A = 40, . . . , 90 meV jzl 
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modification of the crystal field is presently unknown. 
A possible source of the discrepancy between theory and 
experiment is that ARPES is surface sensitive. As shown 
by Pillay et al. 78] , the surface layer of Na ions can give 
rise to an effective t^g crystal field, which substantially 
lowers the e g > bands with respect to the a g bands. 

The results shown in Figure [14] underline the complex 
interplay between kinetic energy and Coulomb repulsion 
in strongly correlated materials. While in some systems, 
such as LaTi0 3 and V 2 3 , orbital polarization tends 
to be enhanced by onsite Coulomb interactions, other 
systems, such as Ca2-;cSr x Ru04 and Na x Co02, exhibit 
more complex patterns. Because of their two-dimensional 
geometry, these oxides comprise subbands of very differ- 
ent widths and spectral shapes. Orbital fluctuations can 
then be enhanced or reduced as a result of Coulomb cor- 
relations. 

As pointed out before, quasi-particle spectra at real 
uj can be derived by analytic continuation of the lat- 
tice Green's function G m (iw„), or by first extrapolat- 
ing Zj rn (iw n ) to real uj and then evaluating G m (u). The 
comparison shown in Figure [T5] demonstrates that both 
methods are consistent. Nonetheless, the latter method is 
seen to retain finer spectral details originating from the 
single-particle Hamiltonian. For instance, the e g spec- 
trum shows two peaks below Ep which evidently are the 
shifted and broadened density of states features near 0.4 
eV and 0.8 eV below the Fermi level. Also, the peak close 
to Ep exhibits some of the fine structure of the single- 
particle density of states. These details are lost if the 
spectrum is instead derived via analytic continuation of 
the lattice Green's function. 



3.3. Degenerate Five-Band Model 

The discovery of superconductivity in iron pnictides 
and chalcogenides has stimulated a variety of investiga- 
tions concerning the role of Coulomb correlations in these 
materials. In contrast to high-T c cuprates, which have 
antifcrromagnetic Mott insulators as parent compounds, 
pnictides are correlated magnetic metals showing signif- 
icant deviations from Fermi-liquid behavior. Moreover, 
as a result of the multi-band nature of pnictides, the in- 
terplay of Coulomb and exchange interactions give rise 
to phenomena not found in cuprates. 

As there is no clear distinction between t2 g and e g sub- 
bands in compounds such as FeAsLaO and FeSe, we have 
recently extended ED DMFT to five bands, where each 
of the impurity orbitals hybridizes with two bath lev- 
els (n s = 15) [34{. Because of the enormous size of the 
resulting Hilbert space, the cluster Green's function is 
evaluated at T — 0, i.e., for f3 — !> oo in Eq. (|8|), while the 
Matsubara grid corresponds to T w 0.01. 

The main result of this study is the identification of 
a spin freezing transition near the nominal n = 6 occu- 




n 



FIG. 16: (Color online) Low-energy scattering rates and spin- 
spin correlations derived within ED DMFT for a degenerate 
five-band model with Bethe lattice density of states (W = 4); 
U = 4, J = 0.75. Solid dots: T = 0.01, empty dots: T = 0.02. 
Upper panel: n s = 15; lower panel n s = 10. 



pancy of the Fe 3d states, i.e., approximately one elec- 
tron away from the Mott phase at half-filling. Towards 
larger occupancies (electron doping), ordinary Fermi liq- 
uid behavior is restored, whereas occupancies lower than 
n w 6 (hole doping) give rise to more pronounced non- 
Fermi-liquid-like deviations caused by a non-vanishing 
scattering rate. Accordingly, the spin-spin correlation 
function C m (r) = {S mz (T)S mz (0)) changes from Pauli- 
like at large doping to Curie- Weiss-like at smaller dop- 
ing. These findings are closely related to the spin-freezing 
transition [47} in the degenerate three-band model near 
ri = 2 occupancy, i.e., about one electron away from the 
Mott phase at half-filling (see Section 3.1). 

Making use of an accurate single-particle Hamiltonian 
i(k), derived from maximally localized Wannier func- 
tions, and constrained RPA interaction parameters cal- 
culated within the same basis [79(, it was also shown 
that, because of larger Coulomb interactions, FeSe lies on 
the non-Fermi-liquid side of the spin freezing transition, 
while FeAsLaO corresponds to a moderately correlated 
Fermi liquid |80j . Deviations from Fermi- liquid behavior 
were recently also shown to occur in the paramagnetic 
phase of a-Fe [§H . 

To elucidate the extent to which these results are re- 
lated to crystal field effects and orbital polarization, we 
consider a fully degenerate five-band model, with semi- 
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FIG. 17: (Color online) Self-energy of degenerate five-band 
model for several occupancies near the spin freezing transi- 
tion; U = 4, T = 0.02. Upper panel: n 3 = 15, two bath levels 
per impurity orbital. Lower panel: n a — 10, one bath level 
per impurity orbital. Note the onsets and low-energy kinks. 



elliptical density of states components of width W = 4 
(W — 4 eV corresponds to typical pnictide compounds). 
Using interaction parameters similar to those in FeSe, 
we find a spin- freezing transition at about n — 6.2 oc- 
cupancy, implying non- Fermi-liquid properties for n = 6. 
The upper panel of Figure [If)] shows the variation of the 
low-energy scattering rate, 7 m = \lmT, m (iuj n —> 0)|, with 
band occupancy for T = 0.01 and T — 0.02, assuming 
two bath levels per impurity orbital. Both Matsubara 
grids indicate a Fermi-liquid to non-Fermi-liquid transi- 
tion near n — 6.2, . . . , 6.4. The spin correlation function 
C m (/3/2) changes from Pauli to Curie- Weiss behavior in 
the same occupancy range. 

Because of the importance of strongly correlated five- 
band transition metals and transition metal oxides, it 
is of great interest to investigate whether ED DMFT 
calculations with only one bath level per impurity or- 
bital can provide useful qualitative information concern- 
ing crucial dynamical properties of the system. Calcula- 
tions for n s = 10 (the dimension of largest sector of the 
Hilbert space is N = 63504) are extremely fast, requiring 
only few minutes per iteration at typical temperatures. 
Thus, for a qualitative exploration of the phase diagram 
they would be highly valuable, before more accurate but 
time-consuming calculations for n s = 15 are carried out. 

The lower panel of Figure [TJI] shows the scattering rate 



and spin correlation function as functions of band oc- 
cupancy for n s = 10, i.e., for only one bath level per 
impurity orbital. Both quantities now exhibit a less 
abrupt change from Fermi-liquid behavior at large n to 
non- Fermi- liquid behavior at small n than for n s = 15. 
Nevertheless, qualitatively, these results also suggest the 
existence of a spin- freezing transition near n — 6.5 occu- 
pancy. Evidently, this transition is not caused by orbital 
polarization. As in the case of the three-band model 
47], it is purely related to many-electron Coulomb and 
exchange interactions at about one electron away from 
half-filling. 

Figure [17] illustrates the influence of the bath size on 
the self-energy. While uncertainties at low frequencies 
are larger for n s = 10 than for n s — 15, the qualitative 
picture is the same in both cases. At large occupancies, 
ImE m (iu n ) ~ uj n for uj n — > 0, consistent with Fermi- 
liquid properties. At about n « 6.5, Im E m begins to ex- 
hibit an onset, indicating the breakdown of Fermi-liquid 
behavior due to the presence of frozen moments. It would 
be very interesting to compare the ED results shown in 
Figs. [16] and [17] with corresponding data obtained via 
continuous-time QMC DMFT. 

The overall consistency of the five-band results for 
n s = 10 with the more accurate ones for n s = 15 might 
be very useful for future applications involving five-band 
materials. The computationally rather straightforward 
ED DMFT approach for n s — 10 could then be used to 
scan a wide range of system parameters, such as tem- 
perature, pressure (band width), crystal field splitting, 
doping concentration, Coulomb and exchange interac- 
tions, etc. For greater accuracy, more demanding cal- 
culations with n s = 15 could then be performed near 
critical values of some of these parameters. One could 
also consider extending finite-temperature ED DMFT to 
/-electron systems with n c = 7 orbitals, using 7 bath 
levels, i.e. n s = 14. 

We close this subsection by pointing out that the self- 
energies shown in Figure [T7] do not only exhibit a low- 
energy onset towards half-filling, but also a kink near 
LJ n « 0.1 eV. Extrapolation of T, m (iuj n ) to real u) reveals 
that this kink is associated with a resonance in Im S m (w) 
which gives rise to a pseudogap in the density of states 
slightly above the Fermi level. In fact, the non-vanishing 
low-energy scattering rate may be viewed as a result of 
this resonance in the self-energy. As shown in Ref. [82| . 
the pseudogap is remarkably stable over a wide range of 
Coulomb and exchange energies, as long as J is not too 
small. The upper panel of Figure [18] illustrates some of 
these quasi-particle spectra, which are derived by ana- 
lytic continuation of the lattice Green's function. More- 
over, the local spin correlation function exhibits a pro- 
nounced maximum at low energies, suggesting that the 
peak in the self-energy corresponds to a collective mode 
induced by spin fluctuations. These features disappear at 
small J (see the spectrum for U = 3 eV, J — 0). This be- 
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FIG. 18: (Color online) Upper panel: Low-energy region of 
spectral distribution of degenerate five-band model (W = 4 
eV) at occupancy n — 6, calculated within ED DMFT (n s — 
10) for several Coulomb energies: (7 = 3 eV (solid red curve), 
U — 4 eV (short-dashed blue curve), and U = 5 eV (long- 
dashed green curve); J = 0.5 eV. Spectra for J — 0.75 eV 
and J — 1 eV are similar. Also shown is the result for U = 
3 eV, J = (dotted magenta curve). The bare density of 
states is indicated by the black dashed curve. Lower panel: 
Fe 3d quasiparticle spectra of FeAsLaO, calculated within ED 
DMFT {n s = 15) for U = 3 eV, J = 0.75 eV H- Solid (red) 
curve: lattice spectrum; long-dashed (green) curve: cluster 
spectrum. The bare 3d density of states is indicated by the 
short-dashed blue curve 17911. 



havior is intimately related to the fact that, at half- filling, 
the system is a Mott insulator for J = 0.5, . . . , 1.0 eV, 
while it becomes metallic for small J. Thus, the Mott 
phase in this multi-band system, and the collective mode 
in ImS m (u) away from half-filling, are a consequence 
of Hund coupling. Qualitatively similar behavior is also 
found in the degenerate three-band Hubbard model dis- 
cussed in Section 3.1. 

The pseudogap induced via Hund coupling found in 
the degenerate five-band model might be relevant also 
for the understanding of the electronic properties of iron- 
based pnictides. As shown in Figure IT51 (lower panel) for 
FeAsLaO, the bare density of states as well as the quasi- 
particle spectrum (obtained via analytic continuation of 
the lattice Green's function) reveal a pronounced mini- 
mum slightly above Ep. (Also shown is the correspond- 
ing cluster spectrum, which agrees well with the lattice 
spectrum in the vicinity of Ep.) Similar features have 
been obtained for a variety of pnictides (see, for instance, 



Ref. The results for the degenerate five-band sug- 

gest that the pseudogap above Ep is a generic feature 
caused by multi-band Coulomb correlations within the 
3d shell and that its existence depends crucially on real- 
istic values of Hund coupling. The paramagnetic quasi- 
particle distribution of actual pnictides should therefore 
consist of a combination of correlation features associ- 
ated with J and signatures related to the bare density of 
states. Evidence for this behavior seems to have been ob- 
served in recent experimental optical data for BaFe2Se2 
30. 



3.4. Triangular Lattice 

The role of spatial fluctuations is of great importance 
for a variety of strongly correlated systems, in particu- 
lar, in the high-T c cuprates, in organic crystals, and in 
graphcnc. In the next three subsections we use the clus- 
ter extension of finite-temperature ED DMFT to eluci- 
date the non-local properties of the self-energy for various 
two-dimensional single-band materials. 

The layered charge transfer salts k-(BEDT-TTF) 2 A:, 
where X denotes inorganic monovalent anions such as 
Cu[N(CN) 2 ]Cl or Cu 2 (CN) 3 , exhibit a remarkably rich 
series of phases as function of hydrostatic pressure, in- 
cluding Fermi-liquid and bad-metallic behavior, super- 
conductivity, as well as paramagnetic and antiferromag- 
netic insulating phases 8(| 87 1. Since these compounds 
have triangular lattice geometries, with equivalent or in- 
equivalent nearest neighbor hopping interactions, a fea- 
ture of special interest is magnetic frustration. Thus, 
long-range magnetic ordering is increasingly frustrated if 
the lattice becomes more isotropic. 

In Ref. [43| a cluster extension of DMFT combined 
with finite-temperature ED was used to study the Mott 
transition in two-dimensional frustrated systems, in par- 
ticular, the variation of the T/U phase diagram with 
the degree of anisotropy. In the isotropic, fully frus- 
trated limit (t' — t), the phase diagram is qualita- 
tively similar to the one obtained within the single-site 
DMFT. On the other hand, for intermediate anisotropy 
t' s» 0.8t, re-entrant behavior was found. Analogous 
QMC CDMFT results [88] also exhibit re-entrant be- 
havior, albeit at higher temperatures. Also, only the 
phase boundary U\ C (T) was determined. More recent 
continuous-time CMC DMFT calculations for t' = 0.44i 
yield T c = 0.13t [89j, i n ap proximate agreement with the 
ED result, T c w O.llt [Hf. 

The unit cells employed in Ref. [43[ consist of three 
sites, and the bath was assumed to contain six levels 
(n s = 9). In Ref. [40] it was shown that, despite this 
small bath size, for T — 0.02 a rather accurate discretiza- 
tion of the bath Green's function is achieved (t = 1, band 
width W = 9). Figure IT51 illustrates how the quality of 
this discretization improves further if three bath levels 
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FIG. 19: (Color online) Discretization of bath Green's func- 
tion components G m ,o(i^n) for isotropic triangular lattice us- 
ing two (n s = 9) and three (n s = 12) bath levels per molecular 
orbital; t = 1, U = 9, T = 0.02. For clarity, the results for 
n s = 9 are shifted up (m = 1) or down (m = 2, 3) by 0.025. 
Upper panel: real part; lower panel: imaginary part. 



per site are included (n s = 12). We consider here only 
the isotropic case t 1 = t, so that in the molecular-orbital 
basis the Green's function is diagonal, as indicated in 
Eq. (fT5")) . The components Go, m (ioj n ) can then be dis- 
cretized independently. Although the fit for n s — 9 is 
already quite good, several more detailed low-energy fea- 
tures are reproduced better for n s = 12. 

Figure [2D] shows the self-energy components for both 
bath sizes. The differences in this case are caused by 
the fact that the Coulomb energy U = 9 is very close 
to the lower boundary U c \ of the coexistence region [43[ . 
Moreover, the value of U& is slightly smaller for n s = 9 
than for n s = 12. Similar differences close to a phase 
boundary were also found in the degenerate three-band 
model discussed in Section 3.1, for instance, in the up- 
per panel of Figure 4. The chemical potential /i = 3 
nearly coincides with the spin freezing transition, so that 
small changes due to different bath sizes can lead to larger 
changes in scattering rates and other dynamical proper- 
ties. Thus, near phase boundaries the size of the bath 
plays a more subtle role. The important point, however, 
is that this merely affects the precise values of certain 
quantities, such as the critical Coulomb energies, and not 
the qualitative aspects of the phase diagram. 

ED DMFT cluster spectra for isotropic and anisotropic 
triangular lattices can be found in Refs. [4(| and [43j j . 



FIG. 20: (Color online) Self-energy components for isotropic 
triangular lattice using two (n s — 9) and three (n s — 12) bath 
levels per molecular orbital; t = 1, U = 9, T = 0.02. Upper 
panel: real part; lower panel: imaginary part. 



respectively. 



3.5. Square Lattice 

One of the most widely studied applications of clus- 
ter DMFT is the two-dimensional Hubbard model for 
the square lattice which is believed to capture some of 
the essential ingredients of the physics of cuprate high- 
temperature superconductors. The parent compounds 
at half-filling are Mott insulators while in the overdoped 
region at large electron or hole concentrations ordinary 
Fermi-liquid behavior is observed. In the underdoped re- 
gion, significant deviations from Fermi-liquid properties 
are found, including pseudogap formation and a strongly 
momentum dependent scattering rate along the Fermi 
surface. Evidently, such effects cannot be described 
within a single-site or local version of DMFT. The most 
natural unit cell for taking into account short-range fluc- 
tuations are 2x2 clusters as shown in the center dia- 
gram of Figure 2. Several groups have applied T = ED 
CDMFT to investigate the non-local properties of the 
single-band Hubbard model for a square lattice 2^- 31 1 . 

Here we focus on the Hamiltonian defined in Eq. (|13| . 
with t = 0.25 and t' = —0.075 as first and second 
neighbor hopping interactions, respectively (band width 
W = 2). In Ref. 4l| we have used finite-temperature 
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ED CDMFT to study various aspects of this model, in 
particular, the collective mode in the (tt, 0) component 
of the self-energy which gives rise to the opening of the 
pseudogap at low doping, the increase of the scattering 
rate with decreasing doping, and the striking asymmetry 
between electron and hole doping caused by the second- 
neighbor hopping t'. For instance, at U = lOt — 2.5 the 
Mott transition driven by electron doping is first-order, 
with characteristic hysteresis behavior of the orbital oc- 
cupancies and double occupancy. In contrast, the transi- 
tion induced by hole doping is continuous. In the latter 
case, upon doping the van Hove singularity at (tt, 0) shifts 
toward the Fermi energy, whereas for electron doping it 
moves away from Ep. 

To implement ED CDMFT for the square lattice, we 
use the diagonal molecular-orbital or plaquette 27J, l91 | 
basis, as indicated in Eq. (fTl)]) . which yields three in- 
dependent Green's function and self-energy components 
(to = 1,...,3). We denote them as T = (0,0), M = 
(tt, tt), and X = (tt,0), where the latter is degenerate 
with Y — (0, tt). Each of the impurity levels hybridizes 
only with its own bath consisting of two levels, as indi- 
cated in the left-hand diagram of Figure 1. The upper 
panels of Figure HH show the discretization of the bath 
Green's function components via the supercluster con- 
sisting of impurity cluster plus bath. The Coulomb en- 
ergy U = 2.5 is chosen so that at half-filling the system 
is a Mott insulator, where U c « 5.6* = 1.4 [90] . The 
chemical potential fi = 0.5 corresponds to about 8 % 
hole doping. At T — 0.01, each component Go iTn (iuj n ) is 
seen to be fitted very well with only five parameters: the 
impurity level, two bath levels, and two hopping terms. 
For completeness, the values are provided in Table I. 

The results shown in Figure [5T] are derived via mini- 
mization of \G a -G$\ as in Eq. (J22), for W n = with 
N = 1. Appendix A. 2 demonstrates that discretizations 
of similar accuracy are achieved for N = and N = 2. 
Also, minimization of \1/Gq — 1/Gq\, as in Eq. (|23|) . for 
W n = l/w^, N = 0, 1, 2, is of nearly the same quality. 

Transforming the molecular-orbital parameters back to 
the non-diagonal site basis, the values given in Table I 
yield the onsite level e = (e\ + £2 + 2e3)/4 = 0.0078, and 
the hopping energies r = — (ei — £2)/^ = 0.2545 ss t, and 
t' = -(ei +e 2 - 2e 3 )/4 = -0.078 « t' . Thus, the impu- 



TABLE I: Impurity levels e m , bath levels ej., and impurity 
bath hopping matrix elements V m ,k for square lattice; n 3 — 
12, U = 2.5, T = 0.01, fj, = 0.5. The corresponding Green's 
function components are shown in Figure [21] 
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FIG. 21: (Color online) Upper two panels: Comparison of di- 
agonal components of lattice bath Green's function Go,m(*w n ) 
(solid and dashed colored curves) and cluster Green's function 
G()\ m (iujn) (symbols) for square lattice, using two bath levels 
per impurity molecular orbital (n s = 12); U = 2.5, T = 0.01, 
/i = 0.5 j4l]]. Lower two panels: Amplified low-energy re- 
gion. Black dotted curves: (?Q m (iu; n ) on a Matsubara grid 
corresponding to T = 0.001. 



rity levels in the diagonal molecular-orbital basis are well 
determined by the asymptotic behavior of Go(ito n ). The 
small deviations from e = 0, r = t, and r' = t' suggest 
that at low u> n a slightly more accurate discretization of 
Go can be achieved than by enforcing these conditions. 

The two lower panels of Figure [5T] show the low- 
energy region of Go(iuj n ) in more detail. The dotted 
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black curves represent the cluster Green's function com- 



1.5 



ponents Gq , t 



on a much finer Matsubara grid, cor- 



responding to T = 0.001. Since cluster spectra are al- 
ways gapped, large oscillations occur at low frequen- 
cies Lo n < O.OI71", in order to accommodate the limit 
ImGg m (iuj n ) — > for w„ — > 0. Evidently, more bath lev- 
els would be required to achieve a satisfactory represen- 
tation of extremely low temperature and low frequency 
properties. [For the hopping 1 ; = 1/4 used here, the tem- 
perature in Fig. 14 of Ref. [30| is T = i/256 = 1/1024.] 
On the other hand, as long as the discussion is limited 
to temperatures of about T = 0.005 or higher, and to 
frequencies larger than about ujq = 7rT, two bath levels 
per impurity orbital are adequate. 

As shown in Ref. 41|, temperatures in the range T = 
0.005, . . . , 0.010 are low enough to provide evidence of a 
collective mode in the (it, 0) component of the self-energy, 
once it is extrapolated from the Matsubara axis to real 
u). This is illustrated in the upper panel of Figure |22| 
which reveals a transition from Fermi liquid properties 
in the overdoped region [6 > 15 %) to pronounced non- 
Fermi-liquid behavior in the underdoped region (S < 15 
%). These self-energies are derived via analytic continu- 
ation from the Matsubara axis to to + ij with 7 = 0.005. 
The collective mode in ImE^ o)^) gives rise to strong 
damping of electronic states near Ep, consistent with a 
pseudogap in the quasi-particle density of states (28. 91|. 
This is also shown in Figure I22| where the middle panel 
represents cluster spectra at hole dopings S — 0.03 and 
S = 0.18 (broadening 7 = 0.02), while the lower panel 
shows the corresponding spectra obtained via analytic 
continuation of the lattice Green's function to real w. 
About 400 . . . 600 Matsubara points are taken into ac- 
count and the same broadening is assumed as in the clus- 
ter spectra. The low-frequency behavior of both spectral 
distributions, in particular, the pseudogap slightly above 
Ep and the concomitant particle-hole asymmetry in the 
underdoped region, is seen to be in good agreement. Far- 
ther below or above Ep, the lattice and cluster spectra 
exhibit large differences. The lattice spectra can also be 
calculated by first extrapolating the self-energy to real 
frequencies and then evaluating Eq. (jT5|) at real oj. The 
results are fully consistent with the spectra derived via 
direct extrapolation of the lattice Green's function. 

Additional cluster spectra as a function of hole doping 
can be found in Ref. 92] , which discusses the correlation- 
induced transfer of spectral weight between low and high 
energies in the two-dimensional Hubbard model. 

A dynamical cluster calculation (DCA) for the square 
lattice was recently carried out for 4 x 4 clusters, using 
Hirsch-Fye QMC as an impurity solver [93| . Figure [23l 
shows the comparison of the (ir, 0) component of the self- 
energy at the lowest Matsubara frequency with the cor- 
responding ED CDMFT results obtained in Ref. [41| for 
2x2 clusters. Both approaches are in good agreement for 
doping 5 > 0.1. The difference at 5 = 0.05 is most likely 
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FIG. 22: (Color online) Upper panel: Low-energy region of 
imaginary part of self-energy Ex=(tt,o) ( tiJ ) obtained via ex- 
trapolation to real frequencies for several hole doping concen- 
trations (broadening 7 = 0.005); U = 2.5, T = 0.005. Middle 
panel: Cluster spectral distributions. Solid red curve: hole 
doping 8 = 0.18 with peak at Ef = 0; dashed blue curve: 
S = 0.03 exhibiting a pseudogap above Ef- Lower panel: 
analogous spectra obtained via analytic continuation of lat- 
tice Green's function 14111. 



related to the singular nature of ImS( ff oj(tu„) close to 
the Mott transition. Of course, the larger 4x4 cluster 
provides a much better momentum resolution, in partic- 
ular, a better distinction between scattering properties 
along the nodal and anti-nodal directions of the Bril- 
louin Zone. Nonetheless, the crossover from Fermi-liquid 
to non-Fermi-liquid behavior is driven by E^o) which is 
seen to be well represented within finite-T ED CDMFT. 
For a systematic study of the variation of the self-energy 
components with cluster size, see Ref. 94] . As also shown 
in Ref. 41| , at lower temperature the Fermi-liquid to non- 
Fermi-liquid transition as a function of doping becomes 
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FIG. 23: (Color online) Imaginary part of self-energy com- 
ponent £x=(ir,o)(* w o) a * first Matsubara frequency as a func- 
tion of doping for U = 1.5, t' = 0. Solid (red) symbols: ED 
CDMFT (n s = 12) results for 2 x 2 clusters at T = 0.01 [HI; 
open (blue) symbols: QMC DCA results for 4x4 clusters at 
T = 0.014 [931. 



sharper, in agreement with the QMC results in [93l. l94j . 

The doping variation of the self-energy shown in Fig- 
ure 1231 with Fermi-liquid behavior in the overdoped re- 
gion (S > 0.2), non- Fermi-liquid properties at lower dop- 
ing (0.05 < S < 0.15), and a Mott phase at half-filling 
(5 = 0), is strikingly similar to the one discussed above 
for the degenerate three-band and five-band models in 
sections 3.1 and 3.3, respectively. In the latter cases, 
Mott insulating phases may occur also at integer occu- 
pancies away from half-filling, but these tend to have con- 
siderably larger critical Coulomb energies 34 , 47 , 8(| . On 
the other hand, at not too large values of U, these sys- 
tems exhibit non-Fermi-liquid properties due to the for- 
mation of local moments, including a finite low-energy 
scattering rate, up to about one electron or hole away 
from half-filling, i.e., \S\ < 0.2,..., 0.3. Only at larger 
doping, ordinary Fermi-liquid behavior is restored. Since 
multi-site correlations in the single-band Hubbard model 
may be formulated within a diagonal molecular-orbital 
basis, where the on-site repulsion is converted into intra- 
orbital and inter-orbital Coulomb and exchange interac- 
tions, the fact that common physical phenomena appear 
in both types of systems is plausible. The parallel treat- 
ment of these materials within DMFT, and the close re- 
lation between key physical properties, suggests that the 
existence of inter-orbital or inter-site interaction channels 
gives rise to fundamentally new physics that is absent in 
the single-band, single-site limit. 



3.6. Honeycomb Lattice 

We have recently used finite-temperature ED in com- 
bination with CDMFT to study the semi-metal to insula- 
tor transition in the honeycomb lattice at half-filling [42| ■ 
The Hubbard Hamiltonian is given by Eq. (1131) , with hop- 
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FIG. 24: (Color online) Comparison of lattice bath Green's 
function Go,m{iun) (solid red curves) and cluster Green's 
function Go jm (iui n ) (dashed blue curves) for honeycomb lat- 
tice using only one bath level per impurity orbital (n s = 12); 
U = 4, T = 0.01. Upper panel: real part; lower panel: imagi- 
nary part. 



ping matrix elements ty = t = 1 for nearest neighbors 
and tij = otherwise. The impurity cluster in which 
spatial fluctuations are accounted for explicitly consists 
of six sites, as illustrated in the right-hand diagram of 
Figure 2. Because of the symmetry of this cluster, the 
lattice Green's function can be diagonalized as indicated 
in Eq. (fT7)l . In the site basis G\\ and G13 = G15 are 
imaginary, while G12 = Gi6 and G14 are real. Thus, the 
diagonal components satisfy G2 = —G\ and G4 = — Gjjj, 
giving two independent complex functions. Analogous 
relations hold for the diagonal components of the cluster 
Green's function and self-energy. Although a T = clus- 
ter calculation with two bath levels per impurity molec- 
ular orbital (n s = 18) is feasible 95[, the aim of this 
section is to demonstrate that, for the honeycomb lat- 
tice at half-filling, adequate results can be obtained even 
if only one bath level per impurity orbital is employed. 
Since for n s = 12 the Hilbert space is still fairly large, the 
level spacing of interacting states is rather dense so that 
a reasonable representation of the low-energy behavior of 
the self-energy can be achieved. 

As discussed in Section 2.4, one criterion for the ac- 
curacy of ED DMFT is the discretization of the lattice 
Green's function Go via the finite bath. In the present 
case, the fc-sum on the right-hand side of Eq. © con- 
sists of only one term, so that each diagonal element 
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FIG. 25: (Color online) Self-energy components Eii(io> n ) 
(i — 1,...,4) of honeycomb lattice obtained within cluster 
DMFT for three Coulomb energies at T = 0.01. For sym- 
metry reasons, En and E13 are imaginary, while E12 and 
E14 are real. Upper t wo p anels: ED [42j; lower two panels: 
continuous-time QMC [96f |. 

Go, m (iuj n ) is fitted with three parameters, the impurity 
level e m , the bath level and the hopping term Vfc, 
where k — m + 6. Figure [M] shows the comparison of 
the lattice and cluster Green's functions for U = 4, i.e., 
close to the Mott transition near U c ~ 3.5. Comparisons 
of similar quality are found at U = 3 and U = 5, i.e., on 
the semi-metallic and insulating sides of the transition. 
Undoubtedly, one reason for the surprisingly good agree- 
ment is the semi-metallic nature at half-filling, which en- 
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FIG. 26: (Color online) Upper panel: low-energy region of 
quasiparticle distribution of honeycomb lattice for several 
Coulomb energies at T = 0.005 42]. The noninteracting den- 
sity of states is indicated by the black dotted curve. Between 
50 and 200 Matsubara points are used to extrapolate the lat- 
tice Green's function to real uj. Lower panel: density of states 
over wider energy range for (7 = 5 and (7 = 7. 

sures that the imaginary part of the lattice Green's func- 
tion vanishes in the limit of small uj n even for U < U c . 
Thus, there is no disparity here of representing a metallic 
lattice Green's function via a finite-size cluster Green's 
function. 

The role of Coulomb correlations in the honeycomb 
lattice has recently also been studied within CDMFT by 
Wu et al. [96] who used continuous-time QMC as im- 
purity solver. Figure US] shows the ED and QMC self- 
energy components in the site basis, where En and £13 
are imaginary, while £12 and £14 are real. The compar- 
ison at three different Coulomb energies ranging from 
the semi-metallic to the Mott insulating phases indi- 
cates that there is very good agreement between the ED 
and continuous-time QMC results. As can be seen in 
Refs. [42| and [96j], similar agreement is found for the 
double occupancies and for the quasi-particle spectra. 

Figure (upper panel) shows the low-energy region 
of the interacting density of states for several Coulomb 
energies. These spectra are derived from an extrapo- 
lation of the local lattice Green's function Gii(iui n ) to 
real w. To illustrate the stability of this extrapolation, 
at each value of U several curves are plotted for 50 to 
200 Matsubara points, with an additional small energy 
broadening of the order of O.lw 2 . (For |w| > 1 the broad- 
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ening is kept constant at 0.1.) At U = 3, a tiny gap or 
pseudogap is seen which is near the limit of what can 
be resolved within ED/DMFT. At U = 4, a full gap 
of width A fx 0.25 has opened. Its width increases ap- 
proximately to A w 0.6 when the Coulomb energy is 
increased to U = 5. This trend is consistent with the one 



found in Rcfs. [96, 97]. The variation of the gap at larger 
U, and the appearance of Hubbard bands, are indicated 
in Figure [26] (lower panel) . These results suggest that 
non-local correlations in the honeycomb lattice induce a 
paramagnetic semi-metal to insulator Mott transition in 
the range U = 3 ... 4, in striking contrast to single-site 
DMFT which yields U c « 10 . . . 13 [H, [99|. 



4. SUMMARY AND OUTLOOK 

The accuracy of DMFT based on finite-temperature 
exact diagonalization has been discussed for a vari- 
ety of strongly correlated materials. The cases in- 
clude three-band and five-band single-site applications of 
DMFT, appropriate to transition metal compounds such 
as Ca 2 Ru0 4 , V 2 3 , LaTi0 4 , Na x Co0 2 , FeAsLaO, etc. 
In addition, short-range correlations in two-dimensional 
single-band Hubbard models are discussed for triangu- 
lar, square, and honeycomb lattices, which represent key 
aspects of organic molecular crystals, high-T c cuprates, 
and graphene. The main focus of this work is on the role 
of the temperature and size of the bath which is used 
to discretize the host lattice surrounding the correlated 
impurity. In contrast to the simpler single-band, single- 
site model, which requires typically 3 to 5 bath levels to 
achieve convergence, the multi-band and multi-site mate- 
rials discussed here are usually well represented using two 
bath levels per orbital or site, as long as the temperature 
is not too low. 

The accuracy of the results is illustrated using three 
different criteria: (i) the convergence of the self-energy 
with bath size, (ii) the quality of the discretization of 
the bath Green's function in terms of a finite cluster, 
and (iii) the comparison with analogous results obtained 
within continuous-time QMC DMFT. For a total of 10 
to 15 impurity and bath levels, the Hilbert space for cor- 
related multi-band or multi-site systems is rather large. 
Thus, because of the indirect coupling among the baths 
pertaining to different impurity levels, the spectrum of 
excited states is very dense, even if only two bath lev- 
els per impurity orbital or site are included. Moreover, 
at temperatures larger than 5 to 10 meV, the bath dis- 
cretization is usually very accurate. At lower temper- 
ature, accordingly more bath levels per orbital or site 
should be employed. For five correlated orbitals within 
local DMFT only one bath level per orbital was found to 
give qualitatively correct results, including the transition 
from a Mott phase at half-filling to non-Fermi-liquid and 
Fermi-liquid behavior for increasing doping. Also, results 



for the self-energy of the single-band Hubbard model on 
the honeycomb lattice at half-filling with only one bath 
level per site were found to be in excellent agreement 
with analogous QMC CDMFT results. 

Since ED DMFT calculations, like those within QMC 
DMFT, are performed along the Matsubara axis, an ex- 
tra step is needed to evaluate the Green's function and 
self-energy to real energies. An attractive feature of ED 
is that the cluster quantities can be calculated directly 
close to the real axis. Moreover, the corresponding lattice 
quantities can be derived via analytic continuation. As 
we have shown for several systems, this extrapolation is 
usually rather stable in the low-energy region, where the 
results agree well with the (broadened) cluster spectra. 

The examples discussed in this paper suggest that 
finite-temperature ED DMFT is a versatile and accurate 
scheme for the description of electronic properties of a 
wide range of strongly correlated materials. It is applica- 
ble at large Coulomb energies and for full Hund exchange. 
Moreover, it is free of sign problems and statistical errors. 
Since it is most suitable at low temperatures ED DMFT 
can be regarded as complementary to continuous-time 
QMC DMFT. Studies at higher temperatures pose no 
problem in principle, but computational times increase 
due to the larger number of excited states that come into 
play. In the future it would be interesting to explore the 
correspondence between ED and QMC at higher temper- 
atures than the ones discussed in the present work, and 
to establish the range of accuracy of both schemes at 
very low temperatures. It would also be very interesting 
to study more d-electron and possibly /-electron mate- 
rials in order to investigate further the accuracy of ED 
DMFT when only one bath level per impurity orbital is 
included. 
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Appendix: Bath Discretization 

One of the key approximations in ED DMFT is the 
discretization of the bath Green's function Go in terms 
of a finite cluster. Whereas Go for the infinite lattice is 
continuous at real uj, the analogous cluster version Gq 
consists of a small number of discrete lines. Nonethe- 
less, at finite temperatures both quantities are smooth 
functions along the imaginary Matsubara axis. To find 
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the bath levels Sk and impurity-bath hopping matrix ele- 
ments V m k , the distance functions in Eq. (|22|) or (l23l) are 
minimized. Evidently, this minimization implies some 
uncertainty due to the different weight functions W n that 
may be used to place more or less weight on the impor- 
tant low-energy region. In this Appendix we discuss in 
more detail the discretization of the bath Green's func- 
tion in the case of the three-band model (see Section 3.1) 
and the square lattice (Section 3.5). 

A.l. Degenerate Three-Band Model 

To illustrate the quality of the bath discretization more 
systematically, we show in Figure [27] the fits obtained for 
the three-band model discussed in Section 3.1 for weight 
functions W n = 1/oj% with iV = 0, 1,2. In all cases, the 
impurity levels are held fixed at e m =i,2.3 = 0, in order to 
focus on the variation of the bath levels. Thus, there are 
four fit parameters per impurity orbital. The two upper 
panels are for the distance function specified in Eq. ([22]) , 
while the lower panels are for minimization according to 
Eq. (|23p . The range of chemical potentials is chosen to 
cover the Fermi-liquid region (/x = 2, n ks 1.5) as well 
as the non-Fermi-liquid region close to the Mott phase 
(ji = 4, n 2.5). The intermediate case (// = 3, n « 2) 
is near the spin-freezing transition. 

Figure [27] demonstrates that all six minimization pro- 
cedures yield fits of good quality, in particular, for N = 1 
and N — 2. Naturally, these fits imply slightly different 
bath parameters £k and V m k ■ As shown in Figure [28j 
the differences between N = 1 and N — 2 are less pro- 
nounced than between N = 1 and N = 0. Also, fit- 
ting Go yields more stable and V m k than fitting 1/Gq. 
Since the bath parameters have only an auxiliary func- 
tion, their variation does not have any physical signif- 
icance. Their sole mathematical purpose is to provide 
an accurate representation of the frequency dependence 
of the lattice bath Green's functions Go,m(*w„), Eq. 
Figure [25] indicates that smaller values of are associ- 
ated with smaller values of V m k- Because of this com- 
pensation, the quality of the various fits in Figure [27] is 
approximately the same. 

The fits shown in Figure [27] are for a fixed set of lattice 
bath Green's functions. We have also carried out com- 
plete self-consistency cycles using the six minimization 
procedures specified above. The resulting self-energies 
are in very good agreement. We therefore conclude that 
the arbitrariness caused by the discretization of the bath 
Green's function is of minor importance in ED DMFT. 
On the whole, iteration to self-consistency is found to 
be most reliable and efficient for the distance function in 
Eq. ([22]), with W n = l/u„. 

A. 2. Square Lattice 
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FIG. 27: (Color online) Discretization of bath Green's func- 
tion for degenerate three-band model as in Section 3.1 using 
two bath levels per impurity orbital (n s — 9). Upper two 
panels: minimization of \Go — Gq 1 \, see Eq. H22[| : lower two 
panels: minimization of |l/Go — 1/G§ [, see Eq. (|23jl . Solid 
red curves: lattice bath Green's function Go,m(iw n ) ( same as 
in Figure 5); remaining curves: cluster Green's functions de- 
rived for weight functions W n = l/w„ with iV = (dashed 
green curves), N = 1 (dashed blue curves), and N = 2 (dotted 
magenta curves). 



In Figure [21] we have illustrated the quality of the 
discretization of the bath Green's function of the 2x2 
square lattice for the three independent molecular or- 
bitals corresponding to V, M, and X. The minimization 
was performed using Eq. (|22p with the weight function 
W n = 1 /w„ . Figure [5S] shows results of analogous dis- 
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FIG. 28: (Color online) Bath levels Sk (top panel) and hy- 
bridization matrix elements V m k (lower panel) for fits in Fig- 
ure I27I Solid dots: fi — 2, empty dots: fi — 3, x symbols: 
fi = 4. Left: minimization of [Go — Gg ! |, see Eq. ([2"2"ll ; right: 
minimization of |1/G - 1/Gg |, see Eq. (J23J). The two bath 
levels are denoted by the red and blue colors. 



cretizations using Eq. (J22J) as well as Eq. ([23]) , for weight 
functions W„ = with A7 = 0,1,2. The impurity 

levels e m are kept fixed at the values determined by the 
asymptotic behavior of Go. m , as indicated in Eq. (TK)1) . 
Thus, there are four fit parameters per molecular orbital 
component Go, m - 

As in Figure [371 there is little difference among the six 
minimizations. Those for N = 1 and N = 2 are slightly 
more accurate than for N = 0. Figure [501 shows the vari- 
ation of the bath levels and hopping matrix elements 
V m k- Although there is considerable variation, in partic- 
ular between N — 1 and N = 0, the net effect on the 
quality of the fit of Go, m using only two bath levels is 
small. Evidently, there is important compensation be- 
tween the magnitudes of Eu and V m k- As in Figure [251 
the differences between the results for N = 1 and N = 2 
are much smaller than for N = 1 and TV = 0. Also, 
fitting Go gives much more stable bath parameters than 
fitting 1 /G . 

In Figs. |2"TI and |2T)1 the bath and cluster Green's func- 
tions are compared since these are the important quan- 
tities used as input in the exact diagonalization. If the 
lattice and cluster hybridization functions [see Eqs. ([24[) 
and (f25]) ] arc compared instead, the differences at large 
oj n are slightly larger since, in contrast to Go, these func- 
tions do not have unit asymptotic weight. 
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FIG. 29: (Color online) Discretization of bath Green's func- 
tion for 2x2 square lattice as in Section 3.5. Solid and 
dashed curves: lattice bath Green's function Go t m{ioj n ) (as 
in Figure [21]) . Upper two panels: real and imaginary parts of 
cluster Green's function &Q >m (iui n ) obtained via minimization 
of | Go - Go\, see Eq. p2j). using two bath levels (n s = 12) 
and weight functions W n = l/tJn witn AT = (+), N = 1 
(o), and N — 2 (x). Lower two panels: analogous fits via 



minimization of |1/G - 1/Gg |, see Eq. J23 
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